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LINEAR  THEORY  OF  BOTTOM  REFLECTIONS 

INTRODUCTION 

The  reflection  of  underwater  explosion  shock  waves  from  the 
bottom  of  the  sea  is  an  interesting,  but  difficult  problem.  The 
mathematical  model  used  today  employs  the  acoustic  approximation, 
namely  the  assumption  of  small  amplitudes  and  a  constant  propagation 
velocity.  This  approximation  permits  the  use  of  the  linear  (acoustic) 
wave  equation  which  has  implicit  solutions.  Certainly,  shock  waves 
from  explosions  are  high  amplitude  (non-linear)  waves,  and  a  strict 
treatment  requires  an  entirely  different  approach. 

The  accuracy  of  the  acoustic  approximation  can,  however,  be 
improved  somewhat  by  replacing  the  sonic  velocity  of  the  media  with 
a  correspondingly  higher  velocity  which  accounts  for  the  increase 
of  the  actual  propagation  velocity  with  pressure.  In  the  same  way 
the  decrease  of  the  pressure  with  distance  can  be  calculated  for  a 
non-linear  rather  than  acoustic  wave.  Although  these  linearizations 
used  with  the  acoustic  method  are  crude  approximations  for  high 
pressures,  much  can  be  learned  about  the  salient  features  of  the 
reflection  process. 

The  calculation  of  the  reflection  of  a  plane,  exponential  wave 
from  a  homogeneous,  plane  bottom  is  relatively  simple  (Arons  and 
Yennie  [l])  .  Since  this  treatment  has  shortcomings,  such  as  the  pre¬ 
diction  of  an  infintely  long  precursor  wave,  solutions  for  spherical  waves 
have  become  of  interest.  The  introduction  of  spherical  waves  considerably 
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complicates  the  problem.  Even  the  highly  simplified  problem  of  the 
reflection  of  a  solitary  spherical,  acoustic  wave  has  a  very  involved 
solution. 

One  of  the  best  treatments  of  the  spherical  wave  is  given  by 
J.  H.  Rosenbaum  [2].  In  his  Naval  Ordnance  Laboratory  report 
Rosenbaum  uses  the  method  of  L.  Cagniard  [3]  to  obtain  integral 
solutions  which  are  suitable  for  numerical  calculations.  Rosenbaum's 
report  is  concise  and  accurate,  but  is  not  very  explicity  written. 

It  leaves  the  reader,  especially  a  beginner  in  the  field,  with  many 
unansvrered  questions.  The  paper  is  also  out  of  print  and  not 
generally  available. 

The  present  paper  has  been  written  for  the  benefit  of  the 
practical  user  who  needs  a  good  understanding  of  the  methods 
involved  in  the  derivation  of  these  equations  to  intelligently  use 
and  manipulate  them  and  be  able  to  attack  more  complicated  problems. 

The  lii.e  of  attack  used  in  this  paper  is  close  to  that  of 
Rosenbaum,  but  many  details  omitted  by  him  are  included  here. 
Rosenbaum's  notation  is  used,  whenever  practical,  to  allow  comparison 
with  his  paper.  Additional  equations  are  derived,  and  a  brief 
discussion  is  given  on  the  use  of  complex  computer  arithmetic  in 
evaluating  the  integral  solution. 
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CHAPTER  I 

SOLUTION  OF  THE  LIQUID  BOTTOM  WAVE  EQUATION 
The  purpose  of  this  Chapter  is  to  develop  from  basic  hydrody¬ 
namic  principles  the  acoustic  wave  equation  in  a  fluid  and  obtain, 
using  the  method  of  Cagniard,  an  integral  solution  for  the  pressure 
from  an  underwater  explosion  over  a  liquid  bottom. 

1-1.  GEOMETRY  AND  HYDRODYNAMIC  RELATIONS 


It  is  assumed  that  a  spherical  acoustic  wave  propagating  through 
a  la^er  of  homogeneous  fluid  (water)  of  thickness  D  interacts  with  a 
bottom  consisting  of  the  lower  half  space  of  either  homogeneous,  non- 
viscous  fluid  or  homogeneous,  isotropic,  and  perfectly  elastic  solid. 
Because  of  the  axial  symmetry,  a  cylindrical-polar  coordinate  system 
is  chosen  with  the  source  of  the  wave  at  the  origin  (r,z,0)  =  (0,0,0). 
Being  rotationally  symmetric,  the  problem  is  independent  of  the  polar 
angle  0.  The  task  is  to  find  the  pressure  at  a  receiver  with  coordi¬ 
nates  (r,z)  as  shown  in  Figure  1-1.1.  The  source  is  at  a  distance  h 
above  the  bottom  and  emits  a  free-water  (direct)  pressure  pulse  given 
by 


Pd  =  0  for  TiR/cj 

Pd  =  F(T  “  R/ci)/R  for  t>r/cx 


(1-1.1) 


The  time  t  is  measured  from  the  instant  of  explosion  t  =  0.  The  slant 
distance  from  source  to  receiver  is  R  =  (rs  +  zz)l/^a ,  and  the  sound 
velocity  in  the  water  is  cx  . 
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The  simplest  case,  that  for  a  liquid  bottom,  will  be  treated 
first.  The  more  complicated  equations  for  a  bottom  with  shear 
strength  (rigid  bottom)  will  be  derived  later  in  Chapter  III. 

We  assume  that  the  water  and  the  liquid  bottom  are  ideal  fluids 
in  the  sense  that  effects  of  viscosity  and  heat  conduction  may  be 
neglected.  We  also  assume  that  the  waves  are  of  small  amplitude, 
that  is,  acoustic  waves. 

The  equation  of  motion  of  a  ideal,  homogeneous  fluid  not  affected 
by  any  external  forces  is  then 

§“  +  (v*v)  v  =  -  £  V  p*  (1-1.2) 

where  P  is  the  excess  pressure,  v  is  the  particle  velocity  vector, 
and  p  is  the  density  of  the  fluid.  For  irrotational  flow  the  particle 
velocity  v  can  be  expressed  byv  =  -  v*<p,  where  cp  is  the  scalar  velocity 
potential.  Equation  (1-1.2)  can  then  be  written 

-  ~  iVv)  +  ~  V  va  =  -  j  VP.  (1-1.3) 

Space  and  time  differentiation  are  interchangeable  and  the 
above  equation  can  therefore  be  integrated  with  respect  to  the  space 
variables.  This  gives  the  Bernoulli  Equation  which  in  the  absence 
of  external  forces  reads 


P  .  1  „a  dt£ 

-  +  j  v  =  dT  ' 


(1-1.4) 
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where  we  have  assumed  that  the  changes  of  p  relative  to  its  initial 
value  p0  are  insignificant.  We  also  have  assumed  that  the  fluid  at 
infinity  is  at  rest  and  has  zero  excess  pressure. 

In  the  acoustic  approximation  the  particle  velocity  v  is  small 
compared  to  the  sound  velocity  c.  Since  the  term  P/p  is  of  the  order 
of  magnitude  vc,  the  term  v*/2  is  small  compared  to  P/p  and  may  be 
omitted  to  give  the  equation  for  the  pressure 


P  * 


(1-1.5) 


The  equation  of  continuity  for  ideal  fluids  is 

§£  +  Pe7*v  -  0  ,  (1-1.6) 

where  we  have  again  assumed  that  the  changes  of  p  relative  to  its 
initial  value  pQ  are  insignificant.  Further  we  neglect  all  dissipative 
processes.  This  means  an  element  of  fluid  cannot  exchange  heat  with 
its  neighbors  and  that  no  energy  dissipation  occurs.  Hence  changes 
in  the  physical  state  of  the  element  take  place  at  constant  entropy, 
and  we  write 


8P 

dT 


(1-1.7) 


where  S  denotes  changes  at  the  constant  entropy  of  the  undisturbed 

fluid.  The  expression  V (dP/dp )g  is  called  the  sound  velocity 

and  is  denoted  by  c.  In  the  acoustic  approximation  changes  of  c 

relative  to  the  initial  value  c  are  ignored. 

0 
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Combining  equations  (1-1.6)  and  (1-1.7)  ,  we  obtain 


-P  V*v 
0 


1 

c0  8  dT 


(1-1.8) 


Substituting  v  and  P  in  terms  of  the  velocity  potential  yields  the 
acoustic  wave  equation: 


V*  <P 


1 

C^  * 


(1-1.9) 


We  now  seek  solutions  of  the  acoustic  wave  equation  which,  with 
the  aid  of  the  Bernoulli  Equation  (equation  (1-1.5)),  can  be  used 
to  construct  the  pressure-time  history  of  a  low  amplitude  underwater 
explosion  shock  wave. 

1-2  WAVE  EQUATIONS  FOR  WATER  AND  BOTTOM,  BOUNDARY  CONDITIONS 

It  has  been  assumed  that  the  motion  in  the  water  and  in  the  bottom 
can  be  represented  by  velocity  potentials  cr  and  ,  the  subscripts  i 
and  a  will  be  used  generally  to  denote,  respectively,  quantities  of 
the  water  and  of  the  bottom.  In  cylindrical  coordinates  cp^  and  cp^ 
satisfy  wave  equations  as  follows: 

in  the  water  |^-  +  £  ff1-  +  55#1-  =■  U-2.1) 

in  the  bottom  fpr8-  +  j  fj*-  +  f^f*-  .  U-2.2) 

where  ci  is  the  sound  velocity  in  the  water  and  c8  is  the  sound 
velocity  in  the  liquid  bottom.  The  particle  velocities  are  then 
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vi  =  -V  /  i  =  1»  2.  (1-2.3) 

To  comply  with  the  physical  picture  of  our  problem,  the 
following  boundary  conditions  must  be  imposed: 

1(a).  At  the  surface  of  the  water,  z  =  -  (D  -  h)  ,  the  excess 
pressure  Pi  must  vanish  since  we  consider  the  water-air  interface 
to  be  a  free  surface. 

2(a)  .  Near  the  source  of  the  explosion  the  effects  of  the 
surface  and  the  bottom  do  not  contribute.  The  pressure  is  then  that 
of  the  direct  wave  given  by  equation  (1-1.1)  as  R  approaches  zero. 

3(a)  .  At  the  bottom-water  interface  pressure  and  normal  velocity 
must  change  continuously  from  one  medium  to  the  other.  This  requires 
Pi  ■  P*  and  dcpj/dz  =  d<?a/dz  for  z  =  h.  Hence  differentiating  with 
respect  to  time,  we  have 


1  a  Pj  _  1_  a  Pa 

p  dz  p  dz  (1-2.4) 

i  a 

4(a)  .  The  disturbance  must  vanish  for  z  -  ®  so  that  P8  approaches 
zero  as  z  goes  to  infinity. 

1-3  LAPLACE  TRANSFORMATION 

To  solve  the  wave  equations  we  use  the  Laplace  transformation 
method.  The  Laplace  transform  f(r,z,s)  with  respect  to  time  t  of 
any  function  f(r,z,T)  is  defined 

00 

f(r,z,s)  =  j'  f(r,z,T)  exp  ( -st)  dT  ,  (1-3.1) 

’  o 
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where  s  is  real  and  is  the  transform  variable  replacing  the  time 

T  . 

Taking  the  transforms  of  equations  (1-2.1)  and  (1-2.2)  and  noting 
that  F(t  -  R/ci )  is  zero  for  t  *  R/cx  ,  we  get 


bP~  r  dr 


0  ,  i  -  1,  2. 


d-3.2) 


Similarly  the  transform  of  the  pressure  «  p ^  ( &tp^/&T )  is 


Pi  =  spi  cp± ,  i  -  1,  2. 


(1-3.3) 


Substitution  of  equation  (1-3.3)  into  equation  (1-3.2)  yields 


0,  i  *  1,  2. 


(1-3.4) 


The  transformed  boundary  conditions  are  as  follows: 

Mb)  At  z  =  -(D  -  h)  ,  the  water  surface,  we  have  P^  *  0. 

2(b)  Near  the  origin  Px  approaches  Pd  =  exp  (-sR/cj)  F(s)/R, 
where  F(s)  denotes  the  Laplace  transform  of  F(t)  . 

3(b)  At  z  =  h,  the  bottom-water  interface,  we  have  Px  »  P*  and 
OPl/^z)/p  =  (dPa/bv)/f>  . 

1  9 

4(b)  As  z  tends  to  infinity  Pt  becomes  zero. 

1-4  SOLUTION  OF  THE  LAPLACE  TRANSFORMED  WAVE  EQUATIONS 

Subject  to  the  boundary  conditions  imposed  on  P^  above,  the 
wave  equation  (1-3.2)  can  be  solved  by  the  method  of  separation  of 
variables.  Assuming  a  solution 
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Pi(r  ,  z,s)  =  Ri(r)  Z^z)  , 


(1-4.1) 


equation  (1-3.2)  can  then  be  written 


1 

r  k. 


(1-4.2) 


At  present  s  is  considered  a  constant.  Since  the  terms  on  the  left 
are  functions  of  r  only  and  those  on  the  right  are  functions  of  z 
only,  we  can  have  equality  only  if  each  is  equal  to  a  constant 
independent  of  r  and  z.  Let  the  constant  be  -X*  .  we  then  have  two 
sets  of  ordinary  differential  eouations: 


rPV  +  +  XsrRi  =0  (1-4.3) 

z V  -  (5la  +  Xa)  Z.  =  0.  (1-4.4) 

l 

Equation  (1-4.3)  is  Bessel  s  differential  equation  of  order  zero. 
It  has  the  solution  =  q^Jo(*r)  which  is  finite  at  r  =  0.  The 
Bessel  function  K0(*r)  is  also  a  solution  of  equation  (1-4.3),  but 
this  function  is  infinite  at  r  =  0.  The  function  J0  is  the  Bessel 
function  of  the  first  kind  and  of  order  zero.  For  our  purposes 
Jo  ( lr  *  Iji  best  defined  by  the  following  integral  representation 
(Dix  l  0)  : 


J0(Xr)  =  Re  —  j  exp  (-iXr  cos  w)  dwj.*  (1-4.5) 


♦Throughout  this  paper  Im  and  Re  denote  respectively  the  imaginary  and 
the  real  parts  of  a  complex  variable  or  function. 
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Equation  (1-4.4)  has  solutions 


fi  eitp[l57>  + 


xa> 


i/a 


;]  +  9i  .xp  [-<£,  +  ^)l/’*]-*Tl-4.6) 


The  constants  ,  g,  ,  and  q^  are  to  be  determined  from  the  initial 
conditions. 

The  "eigenvalue"  X  can  take  on  values  from  zero  to  infinity. 

Since  equation  (1-3.4)  is  linear,  P^(r,z,s)  may  be  written  as  a 
continuous  superposition  of  the  product  in  the  form  of  a 

Fourier -Bessel  integral  as  follows: 

P^  (r ,  z,  s)  *  J  {ai(X)  exp  [(^“a  +  *a)  zj 
o  i 

r  «a  x  /a  -i> 

+  b.  (X)  exp  !  -(-— a  t  Xa)  z  1  f  Jo  (Xr)  XdX  (1-4.7) 

l  i  " 

This  is  a  general  solution  of  the  transformed  wave  equation  for 
cylindrical  symmetry.  To  apply  this  solution  to  the  problem  of  bottom 
reflections,  a^X)  and  b^(X)  must  be  determined  so  that  the  boundary 
conditions  are  satisfied. 

The  above  general  solution  P^  does  not  satisfy  the  boundary 
condition  2(b)  pertaining  to  the  wave  source.  Hence,  the  solution 
is  incomplete  for  the  water  layer .  The  expression  Pi  can  be 
interpreted  as  that  part  of  the  wave  which  results  from  the  bottom 
reflection  but  does  not  include  the  direct  wave.  The  complete 


**Throughout  this  paper  the  symbols  (  )  and  v  denote,  unless 

stated  otherwise,  the  positive  square  root. 
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solution  in  the  water  is  the  sum  of  Px  above  and  the  transform 
of  the  direct  wave  which  satisfies  the  wave  equation  (1-3.4)  and 
the  source  condition  2(b) . 

An  examination  of  the  boundary  condition  4(b)  which  refers  to 
the  signal  strength  at  infinity  shows  that  the  coefficient  aa(X) 
must  be  zero  for  the  wave  in  the  bottom. 

Abbreviating 


R  . 

l 


xa) 


l/s 


(1-4.8) 


the  complete  solution  in  the  water  is 

00 

Px  =  J  {aj^X)  exp(Biz)  +  bx(X)  exp(-Pxz)}  J0(Xr)  XdX  +  ,(1-4.9) 

o 

and  the  bottom  solution  is 

00 

Fa  =  J  {ba  (X)  exp(-Paz)}  Jo(*r)  *<3X.  (1-4.10) 

o 

To  obtain  a  solution  of  the  bottom  reflection  problem  an 
appropriate  expression  for  the  transformed  source  must  be  introduced. 
A  convenient  form  of  the  source  is  that  of  a  spherical  step  wave 
because  by  means  of  the  convolution  integral,  explained  in  Appendix  D, 
solutions  for  any  other  wave  form  can  be  obtained.  The  subscript  A 
will  be  used  when  reference  to  the  step  wave  response  is  made.  We 
choose  the  wave  form  PA  defined  by 

PA  =  R  H(T  "  R/Cl }  '  (1-4.11) 
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where  H(t  -  R/cj )  is  the  Heaviside  unit  step  function  defined 
H (  t  -  R/ci )  =  0  for  t  s  R/cj  and  H (  t  -  R/c* )  ■  1  for  t  >  R/c»  .  The 
expression  R/cj  is  the  arrival  time  of  the  direct  wave.  The  Laplace 
transform  of  the  step  source  is 

PA  =  ^  exp(  -sR/ci )  .  (1-4.12) 

Sommerfeld  has  derived  the  following  identity  for  the  transformed 
step  wave: 


6D 

^  exp  (- sR/cx )  Jo(^r)  exp  (-PxJzj)  3^-  dX  .*(1-4.13) 

0 

This  expression  can  be  combined  with  the  first  integral  of  equation 
(1-4.9).  The  complete  solution  for  the  step  wave  and  its  reflection 
in  the  water  is  then 

00 

Ajtr.z.s)  =  J  jaj  (X)  exp  Ox  z)  +  bx  (X)  exp  (-Px*)}  Jo(*r)*d* 


00 

+  s  J  Jo<Xr>  exP  1*1)  • 


(1-4.14) 


In  the  bottom  the  solution  is 


00 

Aa(r,z,s)  =  J  |ba(X)  exp  (-8az)}  Jo(Xr)  XdX  .  (1-4.15) 


*A  good  derivation  of  this  expansion  is  given  by  Ewing,  Jardetzky, 
and  Press  [5], 
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The  boundary  condition  1(b)  at  the  surface  z  =  -(D  -  h)  requires 
that  Pj  (r,z,s)  =  0  at  z  =  -(D  -  h)  and  hence 

ai  exp  [-Bi  (D  -  h)  "|  +  bi  exp  ^Bi  (D  -  h)  J 

+  exp  [-Bx  (D  -  h)"|  =  0.  (1-4.16) 

At  the  bottom-water  interface  z  =  h  the  boundary  conditions  3(b) 

require  Pj  *=  Pa  and  (dPi/dz)/p  =  (5Pa/3z)/p  and  hence 

l  a 

aj  exp  (Bih)  +  bi  exp  (-Bih)  +  exp  (-Bih)  =  ba  exp  (-Bah) 

(1-4.17) 

and 

[ax  exp  (Bih)  -  bi  exp  (-Bih)  -  exp  (-Bih)^j 

pi  1 

Solution  of  these  three  simultaneous  equations  yields 


(1-4.21) 

where  the  ratio  of  the  density  of  the  water  to  the  density  of  the 


—  B 

=  ^ ba  exp(—  Bah). 

(1-4.18) 
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bottom  is  abbreviated  by 


b  =  p  /p 

x  a 


(1-4.22) 


and  K  is  given  by 


k  =  [(ex/e.)  -  b]/[(ax/3t)  +  bj . 


(1-4.23) 


In  later  discussion  (Appendix  E)  we  will  show  that  K  can  be  interpreted 
as  the  plane  wave  reflection  coefficient. 

The  solutions  Aj  (r,z,s)  and  A^(r,z,s)  can  now  be  written: 

In  the  water 


Jo  (Xr)  X 


-  exp  (-23i  D)  ]  exp  (0iz) 
K  exp  (-201 D) 


exp  [-23i  (D  -  h)  ] [l  -f  K  exp  (-23,  h)  ]  exp  (-Piz) 

1  +  K  exp  (-20* D) 


(1-4.24) 


+  exp  (-Pi  |  z  | )  J-  dX 
and  in  the  bottom 


Jo  eh  J°(xr)x  {(1  +  K)[s*rih<B»  -  l  K  B.)  -  ze,  Dl] 

exp(-03z)}  dX  .  (1-4.25) 


These  are  the  solutions  of  the  transformed  wave  equation  which 
satisfy  the  boundary  conditions  of  our  bottom  reflection  problem  for 


a  step  wave  source 
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1-5.  PEKERIS  RAY  EXPANSION 

Since  we  are  primarily  concerned  with  the  solution  in  the  water 
we  now  restrict  the  discussion  to  A*  only.  After  some  rearranging, 
equation  (1-4.24)  may  be  written 

OD 

Aj  =  J  J0  (Xr)  X  {exp  ( -0i  I  z  j )  -  exp  (8x  z)  (1-5.1) 

0 

+  (exp  (Biz)  -  exp  [-81  (2D  -  2h  +  z)]) 


'1  +  K  exp  (-2Pih)*]], 
.1  +  K  exp  (-28i  D)  Jj 


dX  . 


Note  that  the  term  in  brackets  is  the  only  term  containing  K.  This 
expression  can  be  expanded  in  a  Taylor  series  about  K  exp  (-28XD)  =  0. 
Since  |k  exp  (-28iD)|<l,  this  series  is  uniformly  convergent.  We 
obtain  the  following  expansion  for  Ax  (r,z,s)  in  powers  of  K: 

(1-5.2) 

00 

Al  “  s  J  J°  *i’lX  {k°  £exp  ( — 8x  [  z  | )  —  exp  ( — 8i  [2D  -  2h  +  z])l 

o 

+  K^exp  ( — 8x  C 2h  -  z])  -  exp  (-81  [2D  -  z])  +  exp  (-81  [4D  -  2h  +  z]) 
-  exp  (-81  [2D  +  z])J  +  K*£-  exp  (-8i[2D  +  2h  -  z]) 


+  exp  ( — 8i  [ 4D  -  z]  -  exp  (-81  [6D  -  2h  +  z]) 
+  exp  ( — 8i  [4D  +  z] )  J+  . .  .]•  dX  . 
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Since  this  series  is  uniformly  convergent .  it  can  be  integrated  term 
by  term.  The  above  "ray  expansion"  is  credited  to  Bromwich  [6] 
and  is  used  by  Pekeris  [7].  The  K#  terms  contain  the  direct  wave 
and  the  negative  surface  reflection.  For  each  exponent  n  of  Kn,  n*l, 
there  are  four  terms  which  can  be  represented  by  four  separate  rays, 
denoted  by  m  =  1,  2.  3,  4.  These  rays  will  be  explained  in  more 
detail  below.  Physically,  n  denotes  the  number  of  times  the  ray 
has  been  reflected  from  the  bottom. 

1-6.  IMAGE  SOURCES 

Each  term  in  the  above  series  expansion  can  be  visualized  as  the 

signal  from  an  image  source  located  at  (0,  d  ) ,  where  d  is  the 

n  m  n  m 

height  between  the  gauge  and  the  image  source  (n.  m)  as  shown  in 

Figure  1-6.1.  Each  term  in  equation  (1-5.2)  can  be  written 

(-l)n  +  "V1  exp  (-Pi  ndm)  .  To  find  the  path  of  multiple  reflections 

we  can  make  use  of  the  images  and  the  "images  of  the  images"  with 

respect  to  the  surface  or  the  bottom  as  shown  in  Figure  1-6.1.  We 

must  remember  that  this  method  of  images  is  based  on  the  requirement 

that  the  angle  of  incidence  equals  the  angle  of  reflection. 

From  equation  (1-5.2)  we  can  derive  the  following  formulae  for 

the  four  values  of  d  for  a  given  nils 

n  m  ^ 


d 

=  2(n  -  1)D  +  2h  -  z 

(1-6.1) 

n  i 

d 

n  i 

=  2n  D  -  z 

(1-6.2) 

d 

=  2(n  +  l)D-2h  +  z 

(1-6.3) 

n  3 

d 

n  4 

=  2n  D  +  z  . 

(1-6.4) 
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FIG.  1-6.1  REFLECTED  RAYS  AND  IMAGES 

The  wave  trains  from  the  source  to  the  gauge 
propagate  along  the  rays  shown  as  solid  lines.  To  an 
observer  at  the  gauge  the  trains  seem  to  emerge  from 
the  various  images.  For  instance,  the  wave  train  S-B 
appears  to  originate  at  IMAGE  (1,2). 

The  location  of  IMAGE  (1,2)  can  be  determined  by 
constructing  the  IMAGE  SURFACE  1  which  is  the  water 
surface  mirrored  in  the  bottom  interface.  IMAGE  (1,2) 
is  then  IMAGE  (1,1)  mirrored  in  the  IMAGE  SURFACE  1, 
and  IMAGE  (1,1)  is  in  turn  the  source  mirrored  in  the 
bottom  interface. 
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For  n  =  0  the  surface  reflection  and  the  direct  wave  are  denoted 
respectively  by  m  =  1  and  m  =  2.  The  equations  above  then  reduce 
to  the  following: 

d  =  2D  -  2h  +  z  (1-6.5) 

o  1 

d  =  | z |  .  (1-6.6) 

o  a 

The  significance  of  the  four  values  of  m  becomes  evident  if  we 
use  the  following  notation  (Rosenbaum  [8]) .  Let  B  represent  a 
reflection  from  the  bottom  and  S  represent  a  reflection  from  the 
surface.  Then,  for  n  =  1,  we  have  the  following  reflections  for 
m  =  1,  2,  3,  and  4  respectively:  B,  S  -  B,  S  -  B  -  S,  and  B  -  S. 
Similarly,  for  n  =  2,  there  are  the  following  reflections  for 
m  =  1,  2,  3,  and  4  respectively:  B-S-B,  S-B-S-B,  S-B- 
S  -  B  -  S,  and  B  -  S  -  B  -  S. 

It  is  interesting  to  note  that  for  a  given  value  of  n>l  the  four 
values  of  m  represent  the  combinations  of  alternately  S  and  B 
which  contain  the  letter  B  repeated  n  times. 
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1-7.  INTRODUCTION  OF  AN  INTEGRAL  REPRESENTATION  FOR  THE  BESSEL 
FUNCTION  J0(Xr) 

Each  term  in  equation  (1-5.2)  is  the  Laplace-transformed  step 
wave  pressure  response  resulting  from  image  source  (n,m) .  For 

nal  we  denote  these  transformed  pressures  by  p  .  The  total 

n  m 

transformed  pressure  in  the  water  is  then 

—  _  _  4  •  _ 

Ax (r ,z,b)  -  P '  +  P  +  V  Y  P  ,  (1-7.1) 

oioa  ^nm 

^®1 f 2  0  . . .  n=l f  2  f ... 

where  P  and  P  represent  respectively  the  surface  reflection  and 
0  10  8 

the  direct  wave  terms.  The  general  expression  for  is  then 
obtained  from  equation  (1-5.2) : 

00 

nPm  =  (-l)n+mij  Jo(Xr)  exp  (-5lndm)dX  .  (1-7.2) 

o 

The  Bessel  function  J0(Xr)  may  be  written  as  the  following 
integral: 


tt/2 


Jo(Xr)  *  Re  I  exp  (-iXr  cos  w)  dwl  .  (1-7.3) 

LTT  J  «J 


Substitution  of  this  integral  for  J0(Xr)  in  equation  (1-7.2)  yields 


w  n/2 

"  (-l)n+m  ^  J  ^Re  J  exp(-iXr  cos  w)  dw J 


Pi_1XKn  exp  (“3indm)}dX. 


(1-7.4) 
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Making  the  change  of  variable  *u  ■  X  with  u  real  and  interchanging 
the  order  of  integration,  may  be  written 


P 

n  m 


<-l)n+m  Re 


tt/2 


{II  I 


UOj 


-1 


n 


exp  -  s ( ai  d  +  iur  cos  w)  ^du  dw)  ,  (1-7.5) 

_  n  m  j 


where 


aa 


l  /a 

(u3  +  Cl-3) 

l  /a 

(u3  +  ca-3) 


(1-7.6) 

(1-7.7) 


K  =  (O'!  -  bQ'a)/(or1  +  bo8 )  . 


(1-7.8) 


Although  u  will  be  later  treated  as  a  complex  variable,  u  is  real  in 
the  above  integrals  because  both  s  and  \  are  real.  This  means 
the  integration  above  is  performed  along  the  real  axis  of  the  complex 
u  -  plane.  To  carry  out  the  integration  the  path  will  later  be 
changed  to  pass  over  complex  u. 

1-8.  THE  CHANGE  OF  VARIABLES  (u,t) 

In  equation  (1-7.5)  the  term  in  brackets  suggests  the  change 
in  variables  at  constant  w  defined  by 


t  —  o/\  d  +  iur  cos  w  .  (1-8.1) 

n  m 


The  integral  of  equation  (1-7.5)  then  becomes 
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P  = 
n  m 


(-l)n+m  Re 


K 


n 


(-ST) 


dT^  dw 


} 


(1-8.2) 


In  the  second  integral  the  formal  change  of  variables  (u,t) 

leads  to  an  integration  path  in  the  first  quadrant  of  the  complex 

t  -  plane.  However,  it  is  shown  in  Appendix  C  that  the  path  of 

integration  can  be  shifted  back  to  the  real  axis  of  the  t  -  plane. 

That  is,  the  integration  is  performed  for  real  values  of  t.  This 

variable  t,  defined  by  equation  (1-8.1),  is  shown  in  Appendix  F 

to  be  equivalent  to  the  time  r  when  P  is  evaluated  for  n  =  0, 

n  m 

that  is,  for  the  direct  wave  and  surface  reflection. 

1-9.  THE  INVERSE  LAPLACE  TRANSFORM  OF  nPm 

A  change  of  the  order  of  integration  in  the  above  equation  is 
permissible  and  yields 


P 

n  m 


(-1) n+m  Re 


CD 

£  Jc  - 

Cl 


d 

n  m 


K 


n 


dw  exp  (  —st 


(1-9.1) 


This  integral  may  be  interpreted  as  the  Laplace  transform  of  a  function 
with  the  following  properties: 


n^m 


T  <  Ci 


-l 


d 

n  m 


P  = 
n  m 


(_l)n+m 


uai 


T  i  Ci 


_i 


d  . 
n  m 


(1-9.2) 
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It  is  shown  in  Appendix  B  that  the  press-  e  nPm  is  also  zero 

for  t  <  ,6  where  Ci  _1  £  6  £  c:  -1  R .  The  agnitude  6  is  the 

nm  n  m  n  m  n  m  a  nm 

earliest  arrival  time  of  waves  (n,m)  reaching  the  receiver  and  is 

discussed  in  detail#  along  with  other  characteristic  times  of  the 

pressure  response,  in  Appendix  B.  Therefore,  the  above  expression 

for  ^P^  is  the  desired  solution  of  our  problem.  The  change  of 

variable  (u,t)  and  the  interpretation  of  equation  (1-9.1)  as  a 

Laplace  transform  is  the  vital  step  of  Cagniard's  analysis  which  makes 

trivial  the  inversion  of  the  Laplace  transform  P  and  eliminates 

n  m 

the  integration  to  infinity. 

In  the  following  paragraph  a  further  change  of  variable  is  made 
which  simplifies  the  integral  (1-9.2)  and  its  interpretation. 

1-10.  THE  CHANGE  OF  VARIABLE  (w,u) 

A  change  of  variable  from  w  to  u  at  constant  t  yields  for 


P  =  (-1) n+m  Re 
n  m 


e  C 

-a  r  a  a  a  l/*"| 

where  ux  =  R  -i  Tr  +  d  (t3  -  cx  ”  R  )  i  , 
*  n  m  L  nm  nm  J 


ua  =  d-1  (t3  -  Cl-3  d3  ) 

-  n  m  1  n  m 


i/a 


and  R  =  (r3  +  d  ) 
n  m  n  m 


l/a 


(1-10.1) 

(1-10.2) 

(1-10.3) 


From  equation  (1-8.1)  we  can  write 


(1-10.4) 


^  =  (ir  cos  w  +  — ~^)  / ( iur  sin  w)  =  Uur  sin  w)  j 
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From  equation  (1-8.1)  we  also  obtain 

1  /3  p  -,1  /3 

sin  w  =  (1  -  cos3  w)  =  j^uar3  +  (t  -  oj  ndm^  *  j  /ur  •  (1-10.5) 

Hence  for  t  a  6 


n  m 

(-l)n+m 

Re 

{2_ 

iTTi 

r  -*  Kn 

(1-10.6) 

V 

aT1/3  'l 

[u3r3  + 

(T 

-  d 
n 

m  °>  >  ]  du}  • 

1-11.  ELIMI RATION  OF  THE  OPERATOR  Re 

The  above  expression  for  ^P^  can  be  further  simplified  by 
eliminating  the  operator  Re  which  can  be  accomplished  as  follows: 

Any  complex  quantity  W  has  a  real  part  given  by 

Re  (W)  -  (W  +  ff)/2  ,  (1-11.1) 

where  the  tilde  denotes  the  complex  conjugate.  Let  P  be 

n  m 

denoted  by 


ru» 

P  »Re  f(u)  du. 

n  m  JUj 


We  can  then  write 


n*m 


(1-11.2) 


(1-11.3) 
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where  the  tilde  over  the  integral  denotes  the  complex  conjugate.  In 
general,  the  complex  conjugate  of  an  integral  can  be  written  as 
follows: 


(1-11.4) 


But  for  the  present  case  each  expression  in  f(u)  is  composed  of  sums, 
products,  quotients,  and  powers  of  polynomial  functions  of  u  and  u3 . 
Each  component  function  g(u)  of  f(u)  satisfies  the  relations 


g(u)  =  w  +  iv  (1-11.5) 

g(u)  =  w  -  iv  ,  (1-11.6) 


where  w  and  v  denote  respectively  the  real  and  imaginary  parts  of 
g(u)  .  But  f(u)  also  has  the  constant  factor  1/i  =  -i  which  does  not 
satisfy  equations  (1-11.5)  and  (1-11.6) .  Hence  we  write  f(u)  in 
the  form 


f(u)  =  -i  (W  +  iv)  =  -iW  +  V  ,  (1-11.7) 

where  W  and  V  denote  respectively  the  real  and  imaginary  parts  of 
i  f(u).  The  sum  W  +  iv  satisfies  equations  (1-11.5)  and  (1-11.6)  and 
we  obtain 


f(u)  =  -i(W  -  iv) 


(1-11.8) 
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and 


f(u)  =  -i (W  -iV)  =  iW  +  V. 


(1-11.9) 


Equations  (1-11.8)  and  (1-11.9)  yield 


f(u)  =  -  f(u) 


(1-11.10) 


We  can  then  write 


1  pua  i  pua 

nPm  “  2  J  f(u>  du  -  2  L  f(S)  dTf  * 

Ui  Ux 


(1-11.11) 


Interchanging  the  limits  of  integration  in  the  second  integral  and 
noting  that  u9  »  u,  ,  we  have 


l  i,u*  l 

P  =  4  I  f(u)  du  +  i  I  f(u)  du. 
n  n  2  Ju,  2  Ju, 


(1-11.12) 


These  integrals  can  be  combined  into  one 


,Uj 

P  *  \  J  f(u)  du  , 
n  m  £  Ux 


(1-11.13) 


since  the  integration  variable  u  in  the  second  integral  of  (1-11.12) 
can  be  replaced  by  u.  Substituting  the  expanded  expression  for  f(u) 
we  have 


P 

n  m 


t  <  6 

n  m 


n*m 


|“‘u  o*'1  K"  +  <T  -  ndm  «>>1>r/*du  •  T  1  n6 


-va 


m 

(1-11.14) 
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where  K  *  (or,  -  bor8  )/(<*!  +  b<y«).  Setting 


„  -  [uara  +  (t  -  ndm  or,)3] 


i/a 


(1-11.15) 


the  step  wave  response  P  can  be  written 

n  m 


n^m 


=  0 


t  <  6 

n  m 


n*m 


(-1) 


n+m  1 


1_  f^1  K 
ni  Juiu  or, 


n 


du 


t  i  6 
n  m 


(1-11.16) 


The  above  expression  for  P  is  the  solution  of  the  wave 

n  in 

equation  for  the  reflection  of  a  step  wave 


PA  =  H(t  -  c,"xR)/R  (1-11.17) 

from  a  liquid  bottom.  The  notation  (n,m) ,  introduced  in  paragraph 
1-5#  identifies  the  number  of  the  reflections  from  the  bottom,  n, 
and  a  particular  path  m.  Solutions  for  an  arbitrary  pressure  pulse 
can  be  obtained  from  the  step  wave  solution  by  means  of  the 
convolution  integral  explained  in  Appendix  D. 

The  solution  (1-11.14)  or  (1-11.16)  can  be  evaluated  in  two  ways: 
(a)  by  analytically  separating  out  the  real  part  as  done  by  Rosenbaum 
and  described  in  Chapter  II  or  (b)  by  a  numerical  separation  of  the 
real  part  as  proposed  by  Eichler  and  Rattayya  [9]  and  explained  in 
Chapter  V. 

The  integration  must  be  made  numerically  in  either  case#  and 
the  characteristic  times  explained  in  Appendix  B  must  be  taken  into 
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account.  In  the  first  case  only  real  variables  are  involved  in 
the  final  computation.  The  second  case  requires  a  computer  program 
which  uses  complex  variables  such  as  provided  for  FORTRAN  IV.  The 
first  method  takes  less  computer  time  although  the  machine  program 
is  more  complicated.  On  the  other  hand  it  will  be  seen  in  Appendix  A 
paragraph  A-13  that  the  expressions  for  multiple  reflections  from 
rigid  bottoms  become  so  complicated  that  Eichler  and  Rattayya's 
proposal  becomes  almost  essential. 

1-12  CONCLUDING  REMARKS 

In  the  derivation  of  the  solution  P  ,  Eq.  (1-11.14),  a  number 

n  m 

of  mathematical  details  have  been  bypassed  for  the  sake  of  a  straight¬ 
forward  presentation.  These  details  are  described  in  Appendices  B, 

C,  D...  and  the  reader  interested  in  a  complete  picture  is  urged 
to  read  these  at  this  point. 

Many  of  the  mathematical  operations  use  complex  variables.  For 
the  sake  of  the  reader  unfamiliar  with  the  theory  of  complex 
variables  a  few  salient  points  of  the  theory  are  summarized 
in  Appendix  A. 


28 


NOLTR  69-44 


CHAPTER  II 

THE  LIQUID  BOTTOM  SOLUTION  IN  THE  FORM  OF  ROSENBAUM 
In  this  chapter  we  derive  the  equations  given  by  Rosenbaum  [2] 
for  the  step  wave  pressure  response  when  the  bottom  is  a  fluid. 
General  equations  for  n>l  are  not  derived  in  this  paper.  The  general 
equations  for  a  liquid  bottom  only  are  given  in  Rosenbaum's  paper, 
but  since  they  are  complicated  and  since  corresponding  equations 
for  a  rigid  bottom  would  be  exceedingly  complicated,  it  is  suggested 
that  P  for  higher  ordered  reflections  be  calculated  by  the  method 
explained  in  Chapter  V. 

We  begin  with  equation  (1-11.14)  which  corresponds  to  equation 

a.  X.M 

(15)  of  Rosenbaum's  paper  except  for  the  factor  (-1)  which 

he  neglects  and  a  scale  factor  G  which  he  includes  but  we  set  equal 

to  unity.  To  simplify  the  notation  we  will  usually  drop  the  subscript 

n  (except  in  P  )  ,  and  the  variables  R  ,  d  ,  etc.  will  be  understood 
n  m  mm 

to  mean  R  ,  d  ,  etc.  Also  6  will  be  taken  to  mean  6  .  Even 
n  m  n  m  n  m 

though  it  may  be  confusing,  we  will  continue  to  use  Rosenbaum's 

notation  and  let  t  =  t/R  and  6  =  6/r  . 

m  m  m  m 

2-1.  PRECURSOR  INTEGRAL  (Case  c8  >  c, ) 

When  t  <  ci  ~l  R  the  limits  of  integration  are  imaginary  since 

m  i  V3 

the  square  root  (t  -  Cj  “  R  )  is  imaginary.  The  branch  cuts  for 
the  square  roots  cri  ,  aa  ,  and  y  are  made  respectively  along  the 
imaginary  axis  from  -ici  “x  to  ici  "l  ,  -ic8-1  to  ic8_1,  and  uj  to  tfi 
as  explained  in  Appendix  A.  The  integration  is  performed  on  the 
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Riemann  sheet  on  which  the  radicals  have  positive  real  parts  for 
positive,  real  u.  The  integration  is  then  carried  out  along  the 
imaginary  axis  from  Uj  to  ui  as  explained  in  Appendix  B  and 
diagrammed  in  Figure  2-1.1  below.  It  is  also  shown  in  Appendix  B 


FIG.  2-1.1  INTEGRATION  PATH  FOR  LIQUID  BOTTOM 

PRECURSOR  (c2  >  c]  AND  8<r  <c]Km) 


that  paths  (2)  and  (3)  of  this  figure  from  -ica-1  to  ica-1  make  no 
net  contribution  to  the  integral  since  the  square  roots  <yj  and  aa 
are  real  along  these  paths.  The  pressure  nPm  is  then  the  sum  of 
contributions  from  paths  (1)  and  (4)  which  approach  the  imaginary 
axis  from  the  right  in  the  limit.  The  end  paths  which  join  the 
vertical  paths  (1)  and  (4)  to  the  integration  limits  ui  and  uj  approach 
zero  length,  as  the  vertical  paths  approach  the  imaginary  axis,  and 
do  not  contribute  to  the  integration. 
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Let  P'  denote  the  integrand  of  P  .  Then 
n  m  3  n  m 


<>ui 


-l 


P  «  P'  du  =  (1)  1 

n  in  J  n  ni  „ 

ui 


r-ic8 

P'  du  +  (4)  P  du  . 

*  .  -i  n  m  v  „  n  m 

1C8  Uj 


(2-1.1) 


Since  the  integrand  is  a  function  of  u  only,  we  may  replace  -u  by 
u.  Then  after  interchanging  the  limits  of  integration  on  path  (1) , 


»lCg  I**  — 

P  =  -  (1)  P '  du  +  (4) 

n  m  '  '  J  ~  n  m  J  ~ 

Ui  ux 


-l 


ica 


-l 


P'  du  . 
n  m 


(2-1.2) 


The  labels  indicating  the  paths  are  retained  so  that  the  proper  signs 
can  be  chosen  for  the  square  roots. 

The  square  roots  o-i  and  \  are  real  on  paths  (1)  and  (4)  ,  but 

aa  is  imaginary.  The  imaginary  part  of  ar8  is  positive  on  path  (1) 

and  negative  on  path  (4)  .  Since  aa  is  contained  in  K, 

we  must  determine  the  real  and  imaginary  parts  of  K  to  evaluate  the 

real  part  of  P  . 

n  m 

Since  the  algebra  becomes  increasingly  complicated,  from  this 
point  the  derivation  will  be  specialized  to  the  case  n  =  1.  As  will 
be  explained  in  Chapter  V,  pressure  histories  for  n  >  1  can  be  more 
readily  obtained  using  complex  arithmetic  to  evaluate  equation 
(1-11.14).  It  must  be  remembered  that  the  n  =  1  case  represents 
rays  which  have  been  reflected  once  from  the  bottom  and  that  this 
case,  particularly  for  m  =  1,  makes  the  most  significant  contributions 
to  the  pressure  time  history. 

The  liquid  bottom  reflection  coefficient  K  when  orx  is  real  and 
positive,  but  aa  is  imaginary  and  positive,  can  be  written  using  real 
square  roots  as  follows: 
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K 


-  ib 
+  ib 


+_ 

+ 


Ca 

ca 


1/8 

T7? 


(2-1.3) 


Separating  real  and  imaginary  parts  of  K,  one  obtains 


Im(K) 


Re  (K) 


-2b  «,  [-(U8  +  Ca-3)]1/3 
aj  a  -  ba(ua  +  eg-8) 

tpa  (ua  4  ea-a) 

ot\a  -  ba  (u*  +  Ca 


and 


(2-1.4) 


The  real  part  of  comes  from  Im(K)  .  The  imaginary  part  of 

P  is  identically  zero  since  P  is  the  real  part  of  a  complex 
n  m  n  m 

integral.  (See  equation  (1-10.6) .)  Since  Im(ars)  is  positive  on 
path  (1) ,  but  is  negative  on  path  (4) ,  Im(K)  has  opposite  signs  on 
the  two  paths.  All  other  expressions  in  the  integrand  ^  P^  do  not 
change  in  sign  from  path  (1)  to  path  (4)  .  Hence  p'  also  has  opposite 
signs  on  the  two  paths,  and  the  integrals  of  equation  (2-1.2)  can 
be  combined  into  one  integral.  Equations  (2-1.2)  and  (1-11.16)  then 
yield 

P  =  (_i)m  1  I*  iC*  HJMS)  du  .  (2-1.5) 

l  m  "  J  ~  y 


Making  the  change  of  variable  u  =  ix  and  substituting  Im(K)  from 

above,  P  becomes 
l  m 


-l 


P 

l  m 


(-1) 


m  4b 


f 

Jlm(ui) 


-a,  l/a 


yl«i 


x  dx 
-  c,-a)  T 


(2-1.6) 
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where  y  and  a*  are  expressed  as  functions  of  x  and 


XmlS.)  =  V*  LTr  -  dm(e>'a  V  -  ,S>1/a] 


After  Rosenbaum,  we  introduce  the  following  notation: 


c°s  1  =  =p,  sin  A  =  u)  =  (d"a  -  x3)1/3, 

m  K  m  K 

m  m 

o  *  (cx-a  -  c*-3)1/3,  T  =  5-,  6  =  ~  =  Ca'1  sin  A  +  o  cos  A  , 

m  r  m  R  m  m 

m  m 


M  =  t  cos  A  +  (ci  “a  -  t  a)  x/a  sin  A  , 
mm  m  m  m 


N  =  t  cos  A  -  (ci  “3  -  t  a)l/a  sin  A  ,  and 
mm  m  1  m  '  m 


uu  -  (o  +  M^J/2  +  C(<7  -  M^J/2]  sin  ^ 


Making  the  substitutions  it,  t  ,  o,  sin  A  ,  and  cos  A  in 

m  m  m 

equation  (2-1.6) ,  we  obtain 


jPm  c  (-l)m+1  “2—  J  {tu(oa  -  U)3)1/3  du)|/-[^uua  +  b3(o3  -  id3)  j 


|"(uja  -  cx  "3)  sin3  A  +  t3-2t  u>  cos  A  +  w3  cos3  A  l1/3}  • 
L  *  mmm  m  mJJ 

(2-1.7) 

Substituting  N  ,  M  ,  and  du)  ■  [fT(o  -  M  )/4]  cos  (tt*/2)  df  ,  we  find 
mm  m 


b(o-M  )  *m2  w(o3-ui3)1/3  cos  df 

P«-(-l)m+1  - - 51  J  2 


(2-1.8) 


x  m 


m 


.  [<D^+ba  (a3-u>3)  ] L<d*+N  M  -ui(N  +M  ) 

mm  mm 
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where  we  have  used  N  M  ■  t  3  -  ci  ~3  sin3  l  and  N  +  M  =  2t  cos  l  . 

m  m  m  m  mmin  m 

The  limits  of  integration  must  now  be  determined.  At  x  =  c, “1, 

we  have  w  =  (cj  -a  -  c8-3)1/3  =  c.  Since  id  =  (o  +  M  )/2  +  [(o  -  M  )/2] 

m  m 

sin  rrj/2,  we  get  that  sin  nt/2  =  1  et  u>  ■  o  and  hence  t  .  >  1.  At 

RIZ 

x  m  imtTTj  )  it  is  easily  shown  that  w  =  M^,  and  consequently  that 
*ml  “"1* 

Returning  to  equation  (2-1.8)  ,  it  can  be  shown  using  trignometric 
identities  that 


n* 

cos  -f- 


(uu  -  M^1/3 
Xo“-  ■uu'PTff 


(2-1.9) 


Thus  the  change  of  variable  to  f  was  made  so  that  the  factor  (u>  -  Mm) 
of  \  could  be  eliminated,  removing  the  singularity  of  y  at  u  »  ^ 

(or  u)  «  Mm  and  =  -1) .  Substituting  the  above  expression  for 

cos  tt|/2  and  canceling  (c  -  u>) 1  /a  ( gu  -  M  )*/3  from  numerator  and 

m 

denominator,  we  obtain 


P 

i  m 


(r 


m 


6) 

m 


P 

i  m 


(-1) 


. ,  b(o-M  ) 
m+1  m 


m 


pi  w(o+w)  [l-sin  ]  d(i 

J_1  L  (1-b3)  u)8+aab3]  (uu-N_)  V* 


m 


i 


T 

m 


£ 


ci“l) 


(2-1.10) 


The  above  equation  is  Rosenbaum's  form  of  P  (t  )  for  t  <  cj  ” 

j  m  m  m 

(except  for  the  factor  (-l)m+1  and  the  scale  factor  G) .  When 

t  <  6  ,  we  have  P  (t  )  -  0  since  the  precursor,  the  first  arrival, 
m  in  i  nt  in 

has  not  yet  arrived.  Replacing  the  terms  in  equation  (2-1.10) 
resulting  from  Im(K)  by  corresponding  terms  from  Im(Kn)  one  can 
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obtain  Rosenbaum's  results  for  higher  ordered  reflections  since 
Im(Kn)  has  the  same  sign  dependence  on  as  as  does  Im(K)  . 
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2-2  MAIM  WAVE  INTEGRAL  (Case  c3  >  cx  ) 

The  main  wave  integration  is  quite  similar  to  the  precursor 
integral  except  that  ux  and  ux  are  now  complex.  The  branch  cuts 
for  <*x  and  or*  are  the  same  as  before,  but  the  branch  cut  for  y  is 
made  from  ux  in  the  fourth  quadrant  to  ufx  in  the  first  quadrant. 
The  square  root  y  now  has  a  positive  real  part  left  of  the  cut  and 
a  negative  real  part  right  of  the  cut.  The  initial  integration 
path  is  shown  in  Figure  2-2.1  as  the  line  uxAux.  . 

i  lm(u) 


Re(u) 


FIG.  2-2.1  PATHS  OF  INTEGRATION  IN  THE  u-PLANE  FOR 
T>  C]  _1Rm  (LIQUID  BOTTOM  c2>  c,) 

Since  Re(  y)  is  positive  on  uxAux  ,  but  negative  on  u>Bux  ,  and  since 
orx  and  or*  have  the  same  signs  on  both  paths,  we  can  write 
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P 

n  m 


P'  du 
n  m 


(2-2.1) 


Cauchy's  theorem  allows  us  to  write 


P 


n  m 


1 

2 


f- 


P  ' 
n  m 


du 


P  '  du 
n  m 


(2-2.2) 


As  the  radius  of  &  tends  to  infinity,  P'  has  a  pole  at  infinity 

n  m 

which  in  Appendix  A  is  shown  to  have  the  residue 


R® 


-  (-l)n+m 
in  R 

m 


(2-2.3) 


Since  the  contour  £i  is  in  the  positive  (counterclockwise)  direction, 
we  apply  the  residue  theorem  to  obtain 


P 

n  m 


(-l)n+in 

R 

m 


P' 
n  m 


du  . 


(2-2.4) 


For  more  details  of  the  procedure  used  to  obtain  the  above  equation 
see  Appendix  A. 

Denote  the  integral  on  La  by  I3  .  Then,  dividing  La  into  the  paths 
(1)  ,  (2),  (3),  (4),  (5),  and  (6)  shown  in  Figure  2-2.1,  the  integrals 
along  (5)  and  (6)  between  -ic8-1  and  ica-1  add  to  zero  as  in  the 
precursor  leaving 

*•  -  I1)/!''-:  „Pn  du  +  J1C"  -i  nPm  du  ‘2-2-5> 


+  (3) 


r 


’C8 


-1C! 


-x 

p'  du  +  (4)  T 
-x  n  m  J 


id 


-ic. 


-l 

P  '  du  . 
-l  n  m 


37 


NOLTR  69-44 


Replacing  -u  by  u  on  paths  (3)  and  (4)  and  interchanging  the 
limits  of  integration  on  paths  (1)  and  (4) ,  we  obtain 


Pic8  _1  -ica 

-  (1)  p'  du  +  (2) 

J  .  -i  n  m  .J . 


.ic9 

ici 

.ica 

ici 


ICi 


-x 

p '  du 
-i  n  m 


(2-2.6) 


-l 


,.iCa 

+  (3)  p'  du  -  (4)  P'  du  . 

J  .  -i  n  m  J  .  -l  n  m 

lCi 


-l 


The  labels  indicating  the  original  paths  of  integration  are  retained 

so  that  the  proper  signs  can  be  chosen  for  the  complex  square  roots. 

Again  we  restrict  the  derivation  to  the  case  n  =  1. 

On  paths  (2)  and  (4)  is  real  and  positive,  but  on  paths  (1) 

and  (3)  it  is  real  and  negative.  Let  yx  denote  the  value  of  y  for 

positive  ttj  ,  let  va  denote  the  value  of  y  for  negative  orj  .  As 

in  the  precursor  the  real  part  of  P  comes  from  Im(K) ,  and  the 

n  m 

integrals  on  paths  (1)  and  (3)  can  be  combined  since  Im(K)  ,  given  by 
equation  (2-1.4),  is  positive  on  (1)  and  negative  on  (3)  because  of 
the  imaginary  square  root  cr3  .  Similarly  the  integral  on  paths  (2) 
and  (4)  may  be  added.  On  all  four  paths  and  Va  are  and 

positive.  The  four  integrals  of  equation  (2-2.6)  can  then  be  combined 
into  one  integral  and  Ia  may  be  written 


(-1) 


m 


ic3 


lCj 


-1 


-1 


u[-(ua+c9 

cri[<iiB-b«( 


a .  i  \ 

(u*+Ca"S)  J  Wi  V8  * 


du 


(2-2.7) 


Making  the  substitution  u  =  ix  and  canceling  the  crx's,  we  have 


1.  =  <-um+1 f’.,  ii  -  U  ax  •  (2-2-8) 
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The  following  Rosenbaum  substitutions  are  used:  t  *  t/R  , 

m  m 

ui  *  (Ci”a  -  xa)l/^a,  a  -  (cj  “a  -  c,-a)x/a,  cos  l  -  d  /R  ,  sin  l  ■  r/R  . 

iu  in  in  in  in 

K  =  t  cos  t  ,  and  L  «  (t  8  -  Cx”8)  sin8  l  .  At  the  lower  limit, 
mmmmm*  m 

x  *  cx  l,  uu  becomes  ui  ■  0,  and  at  the  upper  limit,  x  »  c*-1,  uj  becomes 
uu»  =  a.  Thus  the  integral  I8  on  can  be  written 


T  .  ,.m  4b  r°  «(oa  ~ 

r*  “  <-1>  Sr  J  C(l-b^^  +~o« 


m  o 


TbT7 


(2-2.9) 


{cb  -  +  Lmrl/a  - 1(-  +  va  +  Lj-‘/a}  a«. 


Substituting  Ia  above  into  equation  (2-2.4),  we  get  Rosenbaum's 
result  for  t  >  cj  ”x  R  and  ct  >  Cx  : 


m 


P 

x  m 


iriiT1  ^_b 

Rm  1  +  b 


(2-2.10) 


ax  i/a 


.  ,  ,*m+l  2b  I*  wtc-  -  uf)4/-  /r,  xa  .  .  -,-x/a 

+  ‘-1*  Sr-J  1(1  -  b»)5i  ♦  W7  1[(,“  ~  Km  +  hit,3 


m  o 


-  CU  +  Km)a  +  Lm]  l/a}  da:  . 
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2-3  MAIN  WAVE  INTEGRAL  (Case  ca  >  c* ) 

When  the  velocity  of  the  bottom  is  less  than  the  velocity  of  the 
liquid  above  the  bottom,  there  can  be  no  precursor.  The  integration 
paths  are  slightly  different  from  the  case  c*  >  ci  and  are  shown 
in  Figure  2-3.1  below. 


i  Im  (u) 


T>c1'1Rm  (LIQUID  BOTTOM  c}>c2) 

In  comparison  with  the  previous  case,  and  cra  are  real  on  paths 
(5)  and  (6)  ,  and,  as  before,  the  integrals  on  these  paths  cancel. 
Thus  we  can  write 


P 

n  m 


i-Dn*n 

R 


m 


(2-3.1) 
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Let  the  integral  on  path  La  be  denoted  la .  We  then  replace  -u 
by  u  on  paths  (3)  and  (4)  and  interchange  the  limits  of  integration 
on  paths  (2)  and  (3)  to  get 


la 


.ica 

ici 


-l 


ica  -1 

.ica 

P  du  - 

(3)j  t  P  du 

(2 

ici 

J  ici  “x 

(2-3.2) 


+ 


p'du  . 


On  paths  (1),  (2),  (3),  and  (4)  aa  is  real,  but  crx  is  imaginary. 
The  liquid  bottom  reflection  coefficient  for  ox  imaginary  and  positive 
but  cra  real  and  positive  is 


~a \  1 1  /* 


-  bcra 


Im(K)  = 


Re (K)  = 


i[- (ua  +  Cj  ” 

+  bara 

inary  parts 

Of  K, 

we  obta 

2boa [ - (ua 

+  Cl" 

b^oj3  - 

(u3  + 

CX-3) 

~a)  - 

b3aa  3 

and 


bauaa  -  (u* 


(2-3.3) 


(2-3.4) 


Let  yj  and  ya  denote  respectively  the  values  of  y  for  which 

(1)  crx  is  imaginary  and  positive  and  (2)  erx  is  imaginary  and  negative. 
The  square  root  ax  ,  being  imaginary,  is  positive  on  paths  (1)  and 

(2) ,  but  negative  on  paths  (3)  and  (4).  Hence  Re  (K)  is  the  same 
on  all  the  paths,  and  the  contributions  from  Re (K)  cancel. 

In  contrast  to  the  case  c8  >  cx  ,  the  square  root  y  is  now 
complex  so  there  may  be  conti ibutions  from  both  the  real  and  the 
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imaginary  parts.  In  the  contributions  from  Im(K)  the  imaginary  ori's 

cancel  out  the  numerator  and  the  denominator  of  P  so  that  the  real 

i  m 

contribution  to  V  comes  from  Im(y) . 

i  m 

The  integrals  along  paths  (1)  and  (2)  and  along  (3)  and  (4)  may 
be  combined  to  get  an  equation  analogous  to  equation  (2-2.7)  when 
c«  >  ci  •  Since  yx  and  yM  are  complex  conjugates ,  1/yi  and  l/ya  are 
complex  conjugates,  and  there  is  still  symmetry  about  the  real 
axis.  The  ex  '  s  may  be  canceled,  and  we  obtain 

!.  -  (-1)"  -I-(t-) }du  .  (2-3.5) 

The  substitution  u  *  ix  yields 

ca”1 

la  *  (-l)in+1  f  r-T  --a  -  fir  )-Im(i  )}  dx  .(2-3.6) 

Jci“x  Lb  »a  +(x  -Cj  *}]  1  ?*  Ya  J 


The  following  Rosenbaum  substitutions  are  now  used:  t  ■  t/R  , 

mm 


t  8  cos  2  X  +  cj  “8  sina  X  , 
mm  m 


t  a  -  Cj “a  sina  X  . 
m  m 


cos  X  -d  /R  ,  sin  X  «  r/R  ,  D 
m  m  m  in  mm 

E  «  4  (sin8  X  cos8  X  )  t  8(t  8  -  cx  “8)  ,  and  F 
m  m  m'  m  m  *  m 

Also  let  TB  a  (x8  -  Cl-a)l/a  and  o  ■  (ca-a  -  Cj-8)1^8.  The  limits  of 
integration  then  become  u*  a  (cx  ~8  -  ci-8)1^8  «  0  and 
<Sa  *  (ca”8  -  cx_8)x/a  =  <7  .  After  applying  these  substitutions  to 
equation  (2-3.6) ,  we  obtain 


Ia  .  (-!,■*  {iMi  )-Xm(i  )}  d.  .  (2-3.7) 


1/ 


The  term  Im(l/  yx  )  can  be  written 
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1 


[  —  ( uT“  +D)+2t3  cos3/  +2t  iui  cos  t  ] 
in _ m  mm  in 

R  [(1?+D  )a+E  lx/a 
m  mm 


i/s 


(IP+dV+e  i17*'  (2'3'8) 

'  m  w  m 


m  m 


The  complex  square  root  V1^3  can  be  written  V1//a  ■  \v\l^a  {cos[0  (  yi  ) /2] 
+  i  sin  [6  (vi  )/2]},  where 


(-(w*  +  D  )  +  2t  8  cos*  1  > 

6(yi)  =  arccos  { — —3 — - — -3 - 7773 -  } 

1  [ (w8  +  D  )a  +  E  ]1/8  J 


(2-3.9) 


m  m 


Since  sin  (0/2)  *  [(1  -  cos  6)/2]1//a,  we  obtain 

j —  ("Stt3  +D  )+2t8  cos3  t  . 

imv*/*, .i iy«  )a+E  ]1/4[l:cosCarcc°s{'i^'^>a^J^.  }]]lA  . 

R  L  m  mJ  L  2 

(2-3.10) 


Simplifying,  we  obtain 


Im(Vl/a)  »  ^|  {[(«* 
m 


+  D 


a  i1/*  11/8 

»  +ej  +<^-Fm»}  •  <2-3-u> 


The  terms  yi  and  ya  are  complex  conjugates  so  6(ys)  =  0  (Yi )  +  n  and 
Xm(i-)  «  -  lm(i--)  . 

Ya  Yl 


Hence  we  can  write 


Mt*  -hpi-a— 


(2-3.12) 


m  m 


Substituting  equation  (2-3.12)  into  equation  (2-3.7),  we  get 


_  ,  ..m+1  4\£b  r°  w(g*-wa)1/a 

Ia  ”  ttR  J  L  (l-baj  ^+^ab3  J 

m  0 


.=•)*/»  (*>  +Dm) ■  +E  r 1 ' /a  +  ( a*  -Fm)  ,  1  / 


m  m~ 

(aa+D  ) a+E 
x  m  m 


\  '  - 

•J  du>  . 

(2-3.13) 
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Equations  (2-3.13)  and  (2-3.1)  then  give  the  pressure  equation  for 
Ci  >  c»  : 


P  (  t  ) 
i  m'  m 


P  (T  ) 

i  m  m 


(t  <  Cl-) 

m 


(-l)™*1  ,1  -  b 


( 


m 


+  ( 


1  + 

°  l/ 


(  T  >  Cj  _1  ) 

'  m  1 


}m2^b  r  zgtrf  r. 

*  "Rm  J  [(l-ba)i;3+3V;)  «* 


[(^♦Dm)‘+Em]W-Fm>,i 


rrr  m 

(^+D  )a+E 
m  m 


a 

j  d ¥. 


(2-3.14) 


The  above  equation  concludes  the  derivation  of  Rosenbaum's 
expressions  for  Pm  over  a  liquid  bottom.  Similar  equations  for  a 
rigid  bottom  are  derived  in  Chapter  IV. 
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CHAPTER  III 

SOLUTION  OF  THE  RIGID  BOTTOM  WAVE  EQUATION 
The  derivation  of  the  step  wave  responses  Ai(r,z,T)  and 
Aa(r,z,T)  for  a  rigid  bottom  is  quite  similar  to  that  for  a  liquid 
bottom  and  will  be  treated  briefly. 

3-1  WAVE  EQUATION  IN  THE  WATER 

In  the  water  we  choose  a  displacement  potential  0  such  that  the 
displacement  of  an  element  of  fluid  is  s  =  V0.  This  potential  is 
related  to  the  velocity  potential  <p  of  Chapter  I  by  the  equation 
cp  *-30/3t.  Substituting  this  expression  for  cp  in  the  velocity 
potential  wave  equation  (equation  (1-1.9)  and  integrating  the 
resulting  equation  with  respect  to  t,  the  displacement  potential 
satisfies 


i 


Va0 


1 _  d^0 

Cia  9T*  ' 


(3-1.1) 


where 


♦if-* 

r  dr 


al_ 

dza 


(3-1.2) 


and  where  c\  is  the  sound  velocity  in  the  water.  The  potential  0  and 
the  other  potentials  used  in  this  chapter  are  taken  to  be  zero  as  t 
and  R  approach  infinity.  The  radial  and  vertical  displacements  and 
the  pressure  are  then 
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u 


d0 

dr 


H  • and 


.  p  alt , 

P1  dT3  ' 


(3-1.3) 


where  pj  is  the  density  of  the  water. 

3-2  WAVE  EQUATION  IN  THE  RIGID  BOTTOM 

In  the  rigid  bottom*  the  equations  of  motion  can  be  expressed 

in  terms  of  the  displacements  u  and  u  in  the  r  and  z  directions 

r  z 

as  follows  (Ewing,  Jardetzky,  Press  [10]) : 


/dau  .  du  u  dau  \  a*u-\ 

♦  i  5^  -  jf  ♦  -  sssf] 

/ da u  .  du  dau  \  /du  &u  \ 

+  n T  *  -  r  -  SF") 

/»*"«  *\\  „ 

"^dzdr  dra  y  *  Pa  dia 


adta 


(3-2.1) 


(3-2.2) 


where  and  t*  are  Lamp's  elastic  constants  and  p  is  the  density 

a 

of  the  bottom. 

We  can  define  potentials  t  and  U  such  that  the  radial  and 
vertical  displacements  (u  and  u  )  and  the  normal  and  tangential 

IT  Z 

stresses  (T  and  T  )  may  be  expressed  by  the  following  equations: 


u 


r 


d_i  _  du 
dr  ”  dz 


it 

dz 


U 

r 


(3-2.3) 


♦The  rigid  bottom  solution  given  here  is  similar  to  those  given 
by  T.  W.  Spencer  [11],  Eichler  and  Rattayya  [9],  and  Roever,  Vining, 
and  Strick  [12]. 
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zz 


-  +  2*  {0  +  ' 


rz 


“  1*(2SrFz  -  2 


a8u 


_  _JL_  S  Uj 

9za  "r  c4a  a^aj 


(3-2.4) 


(3-2.5) 


where  c3  and  c4  are  respectively  the  velocities  of  propagation  of 
sound  and  shear  waves  in  the  bottom.  Lamp's  constants  X'  and  t*  are 
related  to  c3  and  C4  as  follows: 

p  and  c4  =  (p/P3)l/a  .  (3-2. €} 

Then  the  potentials  +  and  U  satisfy  the  equations 


and 


va*  = 


vau 


u _ i_  a8u 

r7  c4a  dta 


(3-2.7) 

(3-2.8) 


Equation  (3-2.8)  may  be  replaced  by  a  wave  equation  by  using  a 
potential  e;  defined  by  U  =  -a£/dr.  We  then  obtain 


va5 


JL_ 

c4a  3t3 


(3-2.9) 


3-3  BOUNDARY  CONDITIONS  FOR  A  STEF  WAVE  SOURCE 

As  for  the  liquid  bottom  we  now  seek  a  solution  A4(r,z,T)  for 
a  pressure  source  H(t  -  R/cj)/R.  The  boundary  conditions  which  must 
be  satisfied  are  as  follows: 
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1.  The  pressure  Ax(r,z,T)  in  the  water  must  vanish  at  the 
surface  z  =  -  (D  -  h)  . 

2.  The  tangential  component  of  stress  Trz  must  vanisl'  at  the 
interface  z  ■  h. 

3.  The  normal  component  of  displacement  u£  must  be  continuous 
at  the  interface  z  ■  h. 

4.  The  normal  component  of  stress  T  _  must  be  continuous  at 

zz 

the  interface  z  ■  h. 

5.  Near  the  source  the  solution  for  Ax(r,z,T)  must  approach 
PA  -  H(t  -  R/c»)A. 

6.  As  z  tends  to  infinity  t  and  U  must  tend  to  zero. 

3-4  TRANSFORMED  POTENTIALS  J,  J.  U 

Using  the  same  Laplace  transform  method  we  used  in  the  liquid 
bottom  solution  (see  equations  (1-4.14)  and  (1-4.15))*  we  find  that 
the  transformed  solutions  are 

J  {fx  (X)  expOi  z)  +  gx  (X)  sxp(-3xz)}  Jo(*r)  XdX  (3-4.1) 
o 

CD  _ 

-  pTP’J  {J«<Xr>  «xp(-B,  |*|)  £  dX  , 

1  0 

CD 

T  ■  J  g»(X)  exp(-33z)  Jo(^r)  XdX  #  (3-4.2) 

o 

—  QD 

U  -  -  -  J  g4  (X)  exp(-34 z)  |p  [Jo(Xr)]  XdX  ,  (3-4.3) 

o 

where  fx  ,  gx  ,  g*  ,  and  g4  are  functions  of  X  to  be  determined  from  a 
transformed  version  of  the  boundary  data  and  where  3^  ■  (s8c^--  +  X9)1^3 
for  i  ■  1,  3*  4.  The  last  term  of  $  is  obtained  from  the  transformed 
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pressure  source  written  as  a  Bessel  function  integral  (equation 
(1-4.13))  by  using  the  pressure  relation  P  =  -  pda0/&Ta  and  its 
transform  P  =  -  sap^  . 

3-5  TRANSFORMED  STEP  WAVE  RESPONSE  IN  A  RIGID  BOTTOM 

Since  Ax(r,z,s)  is  a  transformed  pressure,  we  can  write  from 
equation  (3-4.1) : 

Aj  =  -sap  0  = 

(3-5.1) 

00  00 

f  |ax  (X)exp(0i  z) +bj  (X)exp(-Sx  z)  }j0  (Xr)  XdX+  J0  ( Xr) exp(-Bx  | z | ) , 

where  ax  (X)  and  bx  (X)  are  respectively  -sap^  f  (X)  and  -8*0^^  (X)  . 

The  application  of  the  transforms  of  boundary  conditions  (1) 
thru  (4)  yields  the  following  equations.  The  requirement  of  the 
vanishing  of  the  transformed  pressure  Ax  (r,z,s)  at  the  surface, 
namely 


Ax  (r,z,s)  =  -  sap  0  =  0  at  z  =  - (D  -  h), 

i 

yields  the  equation 


(1)  fx  (X)exp[-Bx  (D-h)  ]+gx  (X)exp[0x  (D-h)]-  p'  gy^ ■  exp[-Px  (D-h)]»0. 

1  (3-5.2) 

The  vanishing  of  the  tangential  component  of  stress  Trz  and  its 

transform  T  at  the  bottom-water  interface,  namely 
rz  J 


rz 


{= 


M. 


drdz 


-  2 


+  Sl-u  )  = 

dza  c4  a  J 


0  at  z  =  h, 
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yields  the  equation 


2{-9a  Xga  (X)  exp(-9ah)  J0  '(Xr)  }  -  2-[b4  a  Xg4  (  X)  exp(-34h)  J0  '( Xr)  } 


+  | rr  (x94  (X)exp(-04h)  Jo  '(Xr)  ]•  =  0  , 


(3-5.3) 


where  the  prime  in  Jo'(Xr)  above  and  in  Jx  '(Xr),  Ja '(Xr)  ,  and  Jo"(Xr) 
below  denotes  the  derivatives  d/dr  and  da/dra.  Simplifying  the 
above  equation,  we  obtain 

(2)  29s  ga  (X)  exp(-Bah)  +  (2Xa  +  sac4~a)g4(X)  exp(-94h)  =  0.  (3-5.4) 


The  continuity  of  the  normal  component  of  displacement  u  and  its 

z 

transform  uz  at  the  bottom-water  interface,  namely 


“z  =  H  +  I  +  F  “  #  « 


yields  the  equation 


{-Ba  Xga  (X)exp(-9ah)  J0  (Xr)  }  +  {xg4  (X)  exp(-B4h)  J0"  (Xr)  }  + 

{^94  (X)  exp(-94h)  J0 '(Xr)  }  *  {@i  Xfj  (X)  exp(  9i  h)  -  9i  Xgx  (X)  exp(-9jh) 


+  exp(-Bih)}  J0  (Xr)  . 


(3-5.5) 


Simplifying  the  above  equation  using  the  Bessel  function  recurrence 
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formulae  for  J0  '  and  Jj  '  listed  below,  gives 
(3)  -6s  gaU)  exp(-03h)  -  Xag4(X)  exp(-04h) 

(3-5.6) 

=  0ifi(X)  exp (0i  h)  -  0i  gi  (1)  exp(-0i  h)  +  exp(-0ih). 

The  continuity  of  the  normal  component  of  stress  T__  and  its 

“  zz 

transform  T  at  the  bottom-water  interface,  namely 
z  z 

Tzz  =  -  A, (r.z.s)  at  z  =  h, 

leads  to 

xV7  +  2>-{0  +  rtsif®}  -  •*<».*  at  2-h- 

Using  the  transformed  wave  equation  obtained  from  equation  (3-2.7) , 

’*♦  *  if*  T 

we  obtain 

X  '^r\Xg3  (X)exp(-03h)  Jo  (Xr)  ]•  +  2^{c03aXg3  (X)exp(-03h)  J0  Ur)  ] 

+  [-04^4  (X)exp(-04h)  Jo  '(Xr)  -  04  Xg4  (X)exp(-04h)  J0"  (Xr)  ]]•  (3-5.7) 

=  sap^  |x[fi  (X)  exp(0i  h)  +  gi  (X)  exp(-0ih)  ]J0  (Xr)  -  — exp(-0xh)  J0  (Xx)}, 
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Applying  the  recurrence  formulae  for  Jo  '  and  J\  '  below  and  using 
c»  m  C ( A '  +  2^)/Pi3l/#  and  c*  »  (v*/Pt)l^a#  the  above  equation  can 
be  simplified  to  the  following  expression: 

(4)  2pfc**  {|  (2Xa  +  sac*"a)ga(X)  exp(-0,h)  +  Xa04g4(X)  exp(-P4h)} 

«  saPi{fi(X)  exp(0ih)  +  gx  (X)  exp(-Bxh)  -  exp(-pjh)}  . 

(3-5.8) 

Recurrence  Formulae 

Above  we  have  used  the  following  Bessel  function  recurrence 
formulae  [13]: 


and 


Jo  (u)  —  Ja  (u) 

-  2  CJ,(u)], 

(3-5.9) 

Jo  (u)  +  J*  (u) 

»  “  Jj  (u) , 

(3-5.10) 

y.jo  co>]  -  - 

Jl  (u)  . 

(3-5.11) 

Solving  for  fx  (X)  and  gx  (X)  from  (1)  thru  (4)  and  using 
Ax  *  •  we  find 

(X,  -.•pi  f,  (X,  -  £  (3-5.12) 

and 

b.U)  -  g,(X)  -  -eXp[-^aiD  -  h)]|l  :..K  ggffiij}}  .  (3.5.13) 
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where  the  reflection  coefficient  K  is  now  defined 


K  =  9»  [(2Xa+sac4~a)a-4Xagap4]  -  bpas4C4~4 
’  C(2\a+sflC4"a)a-4Xapae4]  +  bes«4C4“4 


(3-5.14) 


3-6  THE  STEP  WAVE  PRESSURE  RESPONSE  P  FOR  A  RIGID  BOTTOM 

n  m 

The  coefficients  aj  (X)  and  bi  (X)  have  exactly  the  same  form 
as  equations  0.-4. 19)  and  (1-4.20)  for  the  liquid  bottom.  The 
function  hx  is  then  in  the  same  form  as  equation  (1-4.24)  .  The 
expression  |k  exp(-23iD) |  can  be  shown  to  be  less  than  one  (Spencer 
[11])  so  that  the  ray  expansion  is  still  valid.  The  derivation 
leading  to  equation  (1-11.14) ,  written  below,  can  then  be  duplicated 
for  a  rigid  bottom  provided  we  use  K  as  defined  by  equation 
(3-5.14)  above  and  make  the  proper  branch  cuts  for  the  multiple-valued 
functions  involved. 

As  for  the  liquid  bottom,  the  step  wave  pressure  response 

P  resulting  from  ray  (n,m)  is 
n  m 


nP  =  °  , 
n  m 


(  t  <  6  ) 

n  nr 


p  .  LOl 
n  m 


K"  Cu*r>  +  (t  -  ndmai,»] 

Ul 


3i-i/a 


du  , 


(3-6.1) 


(t 


a 


6  ) 
n  m 


The  transformation  u  =  X/s  renders  K  in  the  form  necessary 
for  substitution  in  the  above  equation: 


K  =  »i[(2ua+C4~a)a-  4uaoaor4]  -  ba,C4 
ax  [  (2ua+c4  ”a)  a-  4ua<*3  cr4  ]  +  bajc^ 


(3-6.2) 
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where  crj  =  (ua  +  cl~'4)1/^a  ,  aa  =  (u3  +  c3-a)1//a  ,  and  er4  =  (ua  +  c*"3)1^8. 
3-7  STONLEY  POLES 

In  the  following  calculations  we  must  now  consider  the  branch 
points  of  orj  ,  cr3  ,  cr4  and  y  and  the  poles  ±ik  of  the  reflection 
coefficient  K.  The  branch  cuts  are  made  in  the  same  manner  as 
before.  The  poles  ±ik,  called  the  Stonley  poles,  are  the  zeros  of 
the  denominator  of  K.  It  has  been  shown  by  Spencer  [14]  that  the 
Stonley  poles  are  beyond  the  branch  points  ±icj  ,  ±ic3-1  ,  and 
lic^  -1  on  the  imaginary  axis  of  the  u-plane  and  also  that  ±ik  are 
the  only  poles  on  the  top  sheet  of  the  Riemann  plane.  Poles  of 
this  type  indicate  the  existence  of  a  wave.  This  particular  wave 
is  usually  called  a  Stonley  wave  (Cagniard  [15]  refers  to  it  as 
a  Scholte  wave),  and  its  propagation  velocity  is  c  -  1A. 

Dl 

Stonley  waves  propagate  along  the  liquid-solid  interface  in  the 
same  manner  as  Rayleigh  waves  in  a  vacuum-solid  interface.  The 
Stonley  wave  has  been  studied  by  Scholte  [16],  [17]  and  more 
recently  by  Roever,  Vining  and  Strick  [12], 
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CHAPTER  IV 

THE  RIGID  BOTTOM  SOLUTION  IN  THE  FORM  OF  ROSENBAUM 
In  this  chapter  we  follow  the  same  procedure  used  in  Chapter  II 
for  a  liquid  bottom,  and  we  frequently  refer  back  to  equations  in 
that  chapter.  As  before,  we  derive  using  Rosenbaum's  substitutions 
real  integral  solutions  from  equation  (3-6.1),  where  now  K  is  given 
by  equation  (3-6.2).  Since  the  rigid  bottom  K  is  more  complicated 
and  since  it  has  the  Stonley  poles  on  the  top  sheet  of  the  Riemann 
plane,  the  integrals  obtained  will  be  correspondingly  more  complicated 
than  for  the  liquid  bottom. 

4-1  PRECURSOR  INTEGRAL  (Case  c3  >  Ci  >  C4  ) 

The  integration  path  for  <  Ci  ~x ,  which  is  essentially 
identical  to  the  liquid  bottom  path,  is  shown  in  Figure  4-1.1  below. 
(It  must  be  remembered  that  c3  has  replaced  c8  as  the  sound  velocity 
of  the  bottom.)  The  branch  points  ±ic4-1  of  ar4  and  the  Stonley  poles 
±ik  are  outside  the  integration  path. 

For  n  =  1  we  can  now  write  from  equation  (2-1.10) 


i 


P 

m 


(-l)" 


( CTi  -M  ) 

m 


l 


J 


In.(K)  [  1-sin (^  \)  ]  dt 

(«x  -a)l/2(u,-N  )  Vs 
A  n 


(4-1.1) 


where  1  (K)  r.  st  <e  determined  and  oj  -  (cj  -•  c3  *)1/8.  The 
substit  it  ion  of  Im(K)  is  permitted  here  since  the  sign  of  the  integrand 
is  dependent  only  on  the  imaginary  square  root  o3  as  the  liquid  bottom 
solution  was  dependent  only  on  .  The  real  roots  and  ar4  are 
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positive  on  all  the  paths.  After  expanding  the  imaginary  scnaru 
root  a 3,  V  may  be  written  as  follows: 


•o 


K  =  glLi3a!±Sil!iI  -  4iua(-ua-cJ~V/V41  -  bi(-na-ca  “a)1  /*  c4  4 

( 2ua+C4  “a )  3  -  4iua  (-us-ca*a)  1/aa4  1  +  bi(-u--ca  ■*)  i  /8  o4  *4  ' 


(4-1.2) 


where  all  the  square  roots  in  this  equation  are  positive, 
The  imaginary  part  of  K  is  then 


I:  (K)  = 


J>,<2us+c.-Vj[bc.-‘(-ua  -  c.-V/a] 


U,  (2ua+c.  -a)  a  ]a+[  (-ua-c,  _a)  V*  <4uaa,  -be.  •* )  V 


(4-1.3) 


The  substitutions  ox  =  (cj  -a  ~  ca_a)1//a,  oa  =  (|c4  a  -  Ci  ~a  ( ) 1  , 


and  u.  =  (ci -a  -  xa)  ’/a  with  u  *  ix  yield 


Im(K)  = 


-2[.(i 


-a 


-  Cj"3  a. 


^a)aJ^b  (oj8  -  u>a) 


a .  l  /a 


-  c,-a+aa)a  >+{(»,  a-.a)l/a[»(c1-a-.a)(0.a+.a)l/a+  ^]}aj. 


(4-1.4) 


The  above  equations  can  be  further  simplified  by  using 


-a 


A  =  u(%—  -  Ci  "a  +  w3)8, 


-a  ...aw.._3  „a\i/a,.  a  ,  3,i/a 


B  =  u)(cx  a  -  -  a-) 17  3  ( Ga "  +  O 


c  .  f  c«'4 (°ia  -  *a)l/a  . 


Substituting  A,  B ,  and  C  into  Im(K) ,  we  get 


- 
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Im(K) 


-A[b  C4~4(ai8  -  u>V/a 

I  Ml  * 

2 ■  Aa  +  (B  +  C) 


(4-1.5) 


Introducing  the  above  equation  into  equation  (4-1.1)  ,  we  obtain 
Rosenbaum's  result  for  n  =  1  and  Ca  >  Cj  >  c*  s 


P 

1  H’ 


0 


(T 


m 


6  ) 
nr 


P  (  t  ) 

i  rr.  m 


(-Dn+1 


b  (  Ol  -M  )  » 


(oi  +«.) l/*ACl-.in(j  ♦)]  d* 

( u.-Nm)  ‘  / J  [A“  +  (B+C)  *  ] 


(6 


m 


(4-1.6) 


58 


NOLTR  69-44 


4-2  PRECURSOR  INTEGRAL  (Case  ca  >  c*  >  cx ) 

The  case  ca  >  c*  >  cx  is  not  derived  by  Rosenbaum,  but  its 
importance  in  current  work  justifies  treatment  here.  The  integration 
path  (a)  shown  in  Figure  4-2.1  is  chosen  when  |  ux  |  *  C4 ?  and  the 
path  (b)  ,  also  shown  in  Figure  4-2.1,  is  chosen  when  |  ux  |  >  c* -1 . 

The  integral  for  |  ux  |  *  c4  _1  is  almost  identical  to  equation  (4-1.6). 
The  only  difference  is  that  cx -1  >  c* -1  requires  replacing  B  by 
Ba  =  uj(cx*"8  -  u,3)  (ox3  -  w3)l/a(|u,a  -  aa3|)1//a  since  now 
0,  =  dc'3  -  o,-3|)l/3  ■=  (c,-3  -  c'3)l/3. 

Thus  for  n  »  1,  cj  >  c«  >  cx  ,  and  |  ux  |  *  c*’1  we  obtain 


P  (t  ) 
x  m  m 


(t  <  6  ) 
m  m 


m+1 


b(ox  -M  )  S  (ox+u>)l/a  A[ l-sin(x-  *)  ]  d* 


4R  C4 

m 


—  I  \wj  -ru,/  T/  J  «  T 

T-  - r/a-  r  a  - 5*r.  (4-2.1) 

ll  (w-N^1/3  [A3  +  (B,4C)3] 


<a  <  t  <  Cx _I) 

m  m 


When  |  ux  |  >  C4 the  integral  over  paths  (3)  and  (4)  is  again 
zero,  and  ^ is  made  up  of  two  parts:  1.  The  integrals  over  paths 
(2)  and  (5)  whose  sum  is  identical  to  equation  (4-2.1)  except  the  lower 
limit  of  integration  must  be  replaced  by  if  =  f(u  =  ic*-1).  2.  The 

integrals  over  paths  (1)  and  (6)  on  which  era  and  <*4  are  imaginary. 

Derivation  of  the  first  part  of  requires  only  finding  (i  . 

The  variable  ii  is  given  by  the  expression 


W  =  — = -  + 


(cx  -Mj 


m 


sin  tt  t/2  . 


(4-2.2) 


But  at  u  =  ic4 ”1  this  expression  for  w  yields 
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u*  =  (Cl_a  -  c4"a)l/a  =  (|C4-8  -  cj'a|)l/a  =  aa  .  (4-2.3) 


Then  we  obtain 


sin 


”*1 

2 


2oa  -<?!  -M 
_ m 

Oj  -M 
1  m 


and  ii 


2  r 

—  arcsm  l 

TT 


2oa  -Ci  -M 


i_  Ci  -M 


m 


*]• 


(4-2.4) 


To  find  the  second  part  of  P  ,  we  must  examine  K,  the  only 

i  m 

term  that  contains  a3  and  c*4  .  For  or3  and  cr4  imaginary  and  positive, 
denote  the  reflection  coefficient  by  K.  Then  after  expanding  the 
imaginary  square  roots  a3  and  or4  ,  we  obtain 


K  =  («,  [Pua+c,  ~a)a  +  4ua(-us-cJ-,)l/a(-ua-c4-a)‘/a] 

-  ibo1-*(-ua-ca-a),/a)/(«,[(2ua+e4-a)a  +  4ua(-ua-cJ-a),/a(-ua-cr‘)1/4] 
+  ibc*  (-ua-c,  -a) (4-2.5) 

and 

Im(K)  -  (-2bo4"4(-uaC3-a)1//a 

{•,teua+c.-a)a  +  4ua(-ua-ca-a)l/a(-ua-c4-a)1/a]}/)/ 
({«,t(2ua+c,-a>a  +  4ua(-ua-o,-a)>/a(-ua-c.-a)l/a]}a+ 
[bc,'‘(-ua-c,'a)l/a]a/)  .  (4-2.6) 

Make  the  substitutions 

a  =  -  f-i  -s  +  <*>a]3  # 
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B«  -  w(cx  “a  -  u>a)  (ax  a  -  u)*) 1//a  ( |  uja  -  caa|)l/a  , 

and 

C  -  SSLI  ( c xa  -  <«a)l/3  , 


where  uu,  ox  ,  and  og  are  defined  as  above.  Then  divide  the  numerator 
and  denominator  by  16  to  get 

Im(K)=  -  .  (4-2.7) 

2(A-B,)a+Cs 

From  equation  (4-1.2)  we  can  see  that  on  path  (6)  aa  and  cr4  being 
negative  has  the  effect  of  only  changing  the  sign  of  Im(K)  since 
erx  is  positive  on  both  paths  (1)  and  (6).  Hence,  as  before,  the 
integrals  on  paths  (1)  and  (6)  add.  The  second  part  of  can 

m 

thus  be  obtained  from  equation  (2-1.10)  by  substituting  Im(K)  from 
equation  (4-2.7)  above  and  changing  the  limits  of  integration. 

For  n«l,  c»>c«>cj  and  I  ux  |  >  C4  Xi 


,  pm(TJ 

j  m  m 


(t  <  6  ) 
m  m 


P  (T  ) 

1  m  m' 


(-l)m+1b(aj-Mm)  p1  (aj+u))l/aA[l-sin(|  ♦)]  dt 
4  Rm  c44  T»  -Nm)  l/a[A*  +  (Ba+C)*] 

(-l)m+1b(ox-M|n)jt1  (ax  +w) l/a  (A-B» )  [l-sin(^-  ♦)  ]  df 


4  Rm 


(®-Hm)l/"^(A"B»)"  + 


(6  ST  <Ci  X) 
'mm4 


(4-2.8) 
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4-3  MAIN  WAVE  INTEGRAL  (Case  c3  >  Ci  >  C4 ) 

Since  this  case  is  more  involved  than  the  previous  cases  a 

more  detailed  account  will  be  given.  For  t  >  ci  -1  we  have  seen 

m 

**** 

that  uj  and  u*  are  in  the  first  and  fourth  quadrants  respectively, 
and  we  can  replace  the  integral  on  path  ujAui  by 

• 

By  Cauchy's  theorem  we  can  replace  1*  by  I*  ,  I*  ,  Li  ,  and  X«  which 
are  shown  in  Figure  4-3.1.  The  path  Ei  ,  is  again  a  circle  whose 
radius  tends  to  infinity.  Thus,  we  have 

nPm  -  -  1/2  {  f.  X  du  +  h  nPm  du  +  l*„P»  du  +  £  du}'  t4'3'1’ 

Applying  the  Residue  Theorem,  we  obtain 

nPm  =  -ni  [R.  -  <Rlk  +  R.ik)]  -1/2  £»  X  du  ,  (4-3.2) 

where  RM,  R..  ,  and  R  ..  are  the  residues  of  P'  at  ®,  ik,  and  -ik 
ik  -ik  n  m 

respectively  and  are  derived  in  Appendix  A.  The  negative  sign 
immediately  before  +  R_^  appears  because  the  contours  £«  and  2* 
are  taken  in  the  negative  (clockwise)  direction.  The  details  of  the 
above  procedure  are  contained  in  Appendix  A. 

The  branch  cuts  will  be  the  usual  ones.  The  square  roots  <*1  , 
aa  ,  and  n4  have  positive  real  parts  in  the  right  half  plane,  negative 
real  parts  in  the  left  half  plane,  positive  imaginary  parts  in  the 
top  half  plane,  and  negative  imaginary  parts  in  the  bottom  half 
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plane.  The  function  y  has  a  positive  real  part  left  of  the  cut 
from  ui  to  uj  and  a  negative  real  part  right  of  the  cut. 

The  signs  on  the  square  roots  then  require  division  of  the 
branch  line  integral  (again  denoted  by  la)  into  the  ten  paths 
shown  in  Figure  4.3-1.  Since  c3  is  the  largest  velocity,  the  integrals 
between  -ica~x  and  ics -x  cancel.  We  are  then  left  with  integrals 
along  paths  (1)  through  (8) .  Converting  negative  limits  of  integra¬ 
tion  to  positive  and  interchanging  the  limits  of  integration  when 
needed,  la  can  be  written  as  follows: 


ic4 


1C! 


-l 


-l 


p'  du 
n  m 


-  <2> J 


ic4 


-l 


lCi 


-l 


P '  du 
n  m 


-  <3)J 


ica 


ici 


-l 


-l 


P  du 
n  m 


ic3 


-l 


ica 


-x 


ica 


-l 


<4>J  nP*du+(5)J  .  d“  -  <6>J  .np;du 


ici 


-l 


ici 


-x 


ici 


-x 


ic* -x  ic4 

(7>I  nPmdu+<8>J 


icx 


-l 


lCj 


-l 


-x 


P'  du  , 
n  m 


(4-3.3) 


where  the  tags  are  retained  so  that  the  proper  signs  can  be  placed 
on  the  square  roots. 

The  above  paths  may  be  grouped  into  two  sets:  (a)  paths  (1)  , 

(2),  (7),  and  (8)7  and  (b)  paths  (3),  (4),  (5),  and  (6).  In  set 

(a)  the  limits  of  integration  are  ici  “x  to  ic* -x  ,  and  is  imaginary. 

In  set  (b)  the  limits  are  ic3  ~x  to  icj  _x ,  and  ax  is  real. 

Denote  the  reflection  coefficient  in  set  (a)  by  K,  and  let  y^ 
and  7  denote  respectively  y  when  q>i  is  (1)  imaginary  and  positive 
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and  (2)  imaginary  and  negative.  Similarly  in  set  (b)  let  vi  and  ya 
denote  respectively  y  for  <*i  ;  (1)  real  and  positive  and  (2)  real  and 
negative.  Let  la  =  la  +  la' ,  where  I8  is  the  contribution  from  set 
(a)  and  I3  is  the  contribution  from  set  (b) . 

The  integrals  in  set  (b)  correspond  to  the  integrals  of  the 
liquid  bottom  main  wave  (case  ca  >  Ci )  with  c  and  Im(K)  being  replaced 
by  their  rigid  bottom  counterparts  in  equation  (2-2.9).  Equation  (2-2.9) 
can  be  written  for  n  =  1  in  the  form 


(-1) 


m+1 


5SI  I  Im(K){[(»-Km)a+Lm]-l/>-[(U,+Km)’+Lm]-l/,‘}  do,  . 

(4-3.4) 


m 


In  set  (b)  K  is  given  by  equation  (4-1.2)  and  Im(K)  is  given  by 
equation  (4-1.3)  in  terms  of  u  or  by  equation  (4-1.5)  in  terms  of 
uj.  A,  B,  and  C.  The  sign  on  the  integral  is  determined  by  Im(crs)  ; 
hence,  the  integrals  on  paths  (3)  and  (5)  and  on  paths  (4)  and  (6) 
add.  Using  the  substitutions  defined  for  the  precursor,  the 
contribution  to  Is  from  set  (b)  is 


Is 


(-l)m 


b 

ttR  c44 

m 


01 


J 


- 1 _ _ 1 -  ldl„ 

[A3  +  (B+C)a]  \  (w-K  )  a+L  l1/3  [(uH-K  )a+L  I1/3/ 

mm  mm 


(4-3.5) 

Since  a\  is  imaginary,  the  integrals  in  set  (a)  correspond  to 
those  of  the  liquid  bottom  main  wave  (case  cj  >  cs )  with  a  and  Im(K) 
being  replaced  by  their  solid  bottom  counterparts.  Taking  the 
square  roots  crj  ,  as  ,  and  <*4  to  be  positive,  we  can  write  using  only 
real  square  roots 
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K  *  (i (-ua-ci  -a ) l/s  [  (2u3+c*  -a)  a  -  4ius  (-ua-c,  ~a)  l/a0f4  3 

-  bi(-ua-c,-a)l/ac4-*>(i(-ua-ca-a)l/aC(2ua+c4-a)a  (4-3.6) 

-  4iua(-ua-cj‘a)l/aa4:  +  bi(-ua-c3's)l/ac.  ~*)  . 

Then  the  imaginary  part  is 

Im(K)=  (-2[4ua  (-ua-Cl  ~a) l/a  (-ua-c3  ~a)  l/aa4  ]  [b(-ua-Ca  ~a)  l/ac*  "*])/ 

([ 4ua  ( -ua  -Cl  -* ) 1  /a  ( -ua  -c3  "a ) 1  /a a4  ]  a 

+  [(-ua-c1"a)l/a(2ua+c4''a)a+b(-ua-c3"a)l/ac4"4]*/)  .  (4-3.7) 

Making  the  substitutions  u  =  ix,  w  =  (xa-cj -a) 1//a ,  d  =  (ci  ~a-ca  ~a)  l^a , 
and  aa  =  ( |c*  ~a-Cj  -a  | )  x ^a  ,  and  dividing  numerator  and  denominator 
of  equation  (4-3.7)  by  16,  we  obtain 

Im(K)  =  (2[^(cj  _a+u?a)  (o1a+nJ8)x/a(aaa--ar,)l/s][^LIl  (CTl  3+TDa) l/a])/ 

([w(Cl  _a+aa)  ( 0\  a+as)1/a(o3a-wa)  l/a]a  (4-3.8) 

+  -  c  1-a-«»)a4^Il  (ox8+;a)l/aia)  . 

Let  A  =  uj[^ - d"a  -  wa]  , 
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B  =  w(d  "a+u>a)  (ax  2+wS)  l/a(oa8-w3)  l/a  , 

and 

-  bc4  4  /  a  .  -a  x  i  /a 
C  =  — * —  (ox  +ui  ) 


Then  we  can  write 


Im(K) 


_ B  _  Pbc4 

( A  +C )  a  +B  2  L  2 


+ 


Similarly,  the  Re (K)  can  be  written 


(4-3.9) 


Re  (K) 


A^+B8  -C3 
(A+C)  a+B3 


(4-3.10) 


From  equation  (4-3.6)  one  can  see  that  the  sign  on  the  integrand 
^  is  determined  by  or4  and  Im(y)  .  Hence  the  integrals  along  paths 
(1)  ,  (2)  and  along  paths  (7)  and  (8)  add.  The  contribution  Ia 
to  la  can  now  be  written 


If 


(-Dm+l 


c4 


-x 


-1 


(4-3.11) 


Changing  the  variable  of  integration  to  <»,  substituting  Im(K)  from 
equation  (4-3.9),  and  substituting 

[jm(i  ) -  Im(i  )] 

Yi  a 


from  equation  (2-3.12),  the  contribution  Ii'  from  the  set  (a)  is 
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IH 


m+1 

(~l)V2  b 

TTR  C4 
m 


( j^-t-CTi  a)  1</aB 
[(A+S)a+  B3*  ] 


i/a 


{ 


[  (ur*+D  )a+F  1  +(E*-F) 

L  m  n  J  m 

[(^+D  )a  +  El 
m  m 


i/a 

}  du> 


(4-3.12) 


Substituting  and  (R^  +  R  from  Appendix  A  and  I3  and  I» 

from  equations  (4-3.5)  and  (4-3.12)  into  equation  (4-3.2) ,  we  get 
the  pressure  equation  for  n  =  1,  t  >  cj  -1 ,  and  cs  >  Cj  >04: 


P 

1  m 


b(-Dm+1 

2ttR  C4  * 
m 


+ 


VS~  M-l) 

2ttR  C4* 
m 


0 


A(c,a-.a)l/3 

[Aa  +  (B+C)a] 

1— a  .  tv  l/ar 
[  (A+C)a+Bi] 


{ _ 1 _  1 

k  (u>-K  )a+L  ]Va  [  ( ui+K  )  a+L  1^ 


{ 


[  (uf’+D  )a+E  ]l/a  +  (TCa-F  ).i/a 
mm  ml 


[  (uS’a+D  )  a+E  ] 
m  m 


} 


dUJ  , 
(4-3.13) 


where  Ra  =  -  l/rriR  and  the  contribution  L  from  the  Stonley  poles 


with 


(-l)m 


y/2  k 

R  gi 

nr” 


f (aa  +  f)l/a 
*■  aa  +  £ 


(4-3.14) 


ka)a  -  kag. g.  ]  -  [(^  -  **>*  -  *’9.g.] 

ka)  +  2g3g4  +  ka  (2*-  +  2a.)]+  -  -bk~4-  V  (4-3.15) 

9s  9*  4gs  c4  J 


_  _  /va  -tvi/a  _  _  /.a  „  -3vi/a  _  /va  _  -a,i/a 

9i  =  (K  -  Ci  )  ,  g3  =  (k  -  c*  )  ,  g4  =  (K  -  c4  )  , 

a  *  Ta  -  (ka  -  Cj-acosai  ),  and  f  =  4t*  g^cos3; t  . 
m  m  m  3  m 


It  should  be  noted  that  equation  (4-3.13)  differs  from  Rosenbaum's 

results  by  a  factor  of  2  in  each  of  the  two  integrals  of  equation  (4-3.13)  . 
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4-4  MAIN  WAVE  INTEGRAL  (Case  c*  >  c4  >  cx ) 

When  c*  >  c*  >  Cx  the  integration  path  is  that  shown  in 
Figure  4-4.1.  The  branch  cuts  are  the  usual  ones  along  the  imaginary 
axis  and  from  ux  and  .  As  in  the  previous  case,  we  can  write 
using  Cauchy's  theorem  and  the  Residue  theorem 

nPm  -  -ni  tR.  -  <Rik  +  R.ik)]  -jli  .  (4-4.1) 


where 


11  -  I*  nPm  du  • 


The  integrals  on  paths  (9)  and  (10)  again  cancel.  Along  paths 
(1)  ,  (2),  (7),  and  (8)  or3  and  <*4  are  imaginary  and  arx  is  real;  but 
along  paths  (3),  (4),  (5),  and  (6)  crx  and  cr4  are  real  and  only  cr* 
is  imaginary.  Using  the  procedure  and  notation  used  in  the  previous 
case  (c*  >  Cx  >  c4 )  the  integrals  along  paths  (1)  through  (8)  can 
be  reduced  to  the  following  for  n  «  1: 


2(-l)“+1 

l*ic4  ~l 

Im(K)  udu 

r^-  - 

rr 

J  icx  _1 

0*1 

Lyi 

2(-l)m+1 

pic*  "x 

Im(K) udu 

ri.  - 

TT 

J  2  _  -1 

1C4 

L\ 

-y 

-H 


(4-4.2) 


where  K  is  given  by  equation  (4-1.2)  and  K  is  given  by  equation 
(4-2.5)  .  Since  c4  >  Cx  we  must  replace  B  in  K  by  Be  . 
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The  second  integral  above  is  given  by  equation  (4-3.5)  except 
the  lower  limit  of  integration  must  be  changed  to  ca  =  ( |c4  ~a-ci  -a  | )  , 

the  value  of  m  at  u  =  ic4 -1 ,  and  B  must  be  replaced  by  Ba  . 

The  first  integral  of  Ia  is  also  quite  similar  to  equation  (4-3.5). 
The  upper  limit  of  integration  must  be  changed  to  oa ,  and  the  imaginary 
part  of  the  reflection  coefficient  must  be  Im(K)  which  is  given  by 
equation  (4-2.7). 

Since  the  residues  are  unchanged,  we  may  then  write  for 
>  ci  -1  and  ca  >  c4  >  ci 


P  *  5-  +  A 
i  m  Rm  i  m 


(-l)m+1b 

r°* 

(oi  a-(«a )  1  /a  (A-B.  ) 

2nR  c* 4 
n\ 

J 

0 

[  (A-Ba  )  3  +  Ca] 

<-ir+ib 

f01 

A(a,a-.a)‘/a  f 

2nRm  c44 
m 

'‘o. 

[A*  +  (B8  +C)  a  1  '“[(w 

m  m‘ 


-m 


i - ) 


dai 


m'  m' 


- 1 - du<  , 

C(w+K)*+L  ]l/aJ 


m  m 


(4-4.3) 


where  A  is  given  by  equations  (4-3.14)  and  (4-3.15) . 
i  m 

The  above  equation  completes  the  derivation  of  expressions  for 
P  in  Rosenbaum's  form  for  the  most  important  cases  of  a  rigid 

i  m 

bottom.  An  alternative  method  of  evaluating  P  which  is  applicable 

n  m 

for  all  values  of  n  is  discussed  in  the  next  chapter. 
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CHAPTER  V 

THE  USE  OF  COMPLEX  ARITHMETIC  TO  CALCULATE 
PRESSURE  TIME  HISTORIES 

Although  Rosenbaum's  equations  are  in  a  form  quite  suitable 
for  numerical  calculations,  some  difficulty  is  encountered  when 
one  wishes  to  extend  the  solid  bottom  solution  to  n  >  1  because  of 
the  complexity  of  the  Stonley  residues.  Also,  for  all  cases  n  >  1 
most  of  the  advantage  gained  by  separating  real  and  imaginary  parts 
is  already  lost  since  the  equations  are  so  complicated.  We  then 
consider  evaluating  equation  (1-11.14)  using  complex  computer 
arithmetic,  as  suggested  by  Eichler  and  Rattayya  [9]. 

5-1  INTEGRATION  PATHS  IN  THE  COMPLEX  u-PLANE 
Equation  (1-11.16) , 


P 

n  m 


,  . .  n+m  .ux 
«i 


du 


(5-1.1) 


where  y  =  [uara  +  (t  -  d  ori)8]1^8  and  ux  ■  R  -a[-iTr  +  d  (t8  - 

n  m  nm  nm 

cx  ~anRm  can  be  evaluated  using  complex  computer  arithmetic 

if  the  integration  is  divided  into  the  two  time  periods  t  <  cx  1nRm 

and  t  >  cx  -1nRm  which  result  because  of  the  square  root  (t8-cx  anR^) 1//# 

in  ux .  In  either  case  we  can  evaluate  the  integral  over  half  the  range, 

that  is  from  0  to  Ui,  and  double  the  result  to  get  P  since  the 

n  m 

integrand  is  symmetrical  about  the  real  axis.  We  then  integrate 
along  paths  on  which  the  real  part  x  of  u  is  constant,  nd  the 
imaginary  part  y  varies  from  yx  to  y8 .  That  is,  we  integrate  from 
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u  ■»  x  +  iyi  to  u  =  x  +  iy3  ,  and  the  expression  for  is  as 
follows: 


P 

n  m 


2 ( -1) n+r 
rri 


x+iya 

u  Qfx  1  Kn>  1  • 

x+iyi 


The  limits  of  integration  are  as  follows: 


For  t  <  ci  R  ,  x  =  0,  =  cs  ,  and 

*  n  m  * 


R“a 
n  m 


n  m 


For  t  >  x  =  Re(uj)  =  nR 


-a 


R~a 

n  m 


-a 

Ra  - 

t')'/ 

n  m 

d 

n  m 

Cl 

1 

o 

1 

N 

Tr . 

n  m 


(5-1.2) 


5-2  COMPENSATED  INTEGRAND 

When  evaluating  the  integral  above  numerically,  care  must  be 
taken  to  avoid  the  singularity  at  u  =  Uj .  One  of  the  best  ways  of 
handling  such  a  problem  is  the  "compensated  integrand"  which  we 
illustrate  with  the  example  below. 

Denote  the  value  of  K  at  u  =  u:  by  Ki  .  We  can  then  write 


P 

n  m 


IX  +  la 


+ 


2(-i)n+m 

TT  i 

2(-l)n4Tn 

ni 


x+iya 

j  u«i  ”l  (Kn  -  Kjn)  v“l 
x+iyi 

x+iya 

Kj  n  J  UOTi  _1  v"1  du  . 
x+iyj 


du 


(5-2.1) 


The  integrand  of  Ix  is  then  finite  at  all  points  of  the 
integration  path  for  t  ^  ci_1nRm  8ince  the  limit  as  u  -  ux  is  zero. 
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This  integral  can  be  reduced  to  integration  over  a  real  variable 
if  we  denote  the  integrand  P;  =  a  +  ib  and  the  integration  variable 
u  =  x  +  iy.  Then  lx  can  be  written  in  real  integrals  as  follows: 


The  first  and  fourth  integrals  are  zero  since  we  are  integrating 
along  a  path  of  constant  real  part  x  of  u.  The  real  contribution 
to  Ij  is  then  the  second  integral,  and  we  can  write 

X,  =  f  in,  (K»  -  Ov-}  ay  •  (5-2.3) 

Jyi  1  J 

Eliminating  the  "i",  Ix  becomes 


I»  =  2l-lln'Hn  fYa  Re  {uoTj-1  (Kxn  -  Ki  Hjv"1}  dy  .  (5-2.4) 

"Yi 

The  integral  Ia  can  be  integrated  analytically  as  shown  below: 
The  substitution  iu  =  (ex  _a  +  ua)1//a  in  I3  of  equation  (5-2.1) 
yields 


x.  = 

*  TTl 


n+m 


n  r 

K.  J 

«> 1 


r  *• 

Ln  m 


O) 


-2  d 
n  m 


T  U)  + 


(TS-C,-2 


dou,  (5  2.5) 


where  ti*  =  [cj  “a  +  (x  +  iyx)a1x//a  and  u*  =  [ex  ”a  +  (x  +  iya)a]1//a. 

The  expression  in  brackets  is  of  the  form  1/y/x  where  X  =*  awa  +  bu>  +  c. 
The  integrand  is  then  the  derivative  of  the  function 


F(u,)  =  Ja  log  +  aVT  +  ~  }, 


(5-2.6) 
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or  in  expanded  form, 


F(uO  =  log  {[  R3  a3 -2  d  tu)  +  (t3-^  ~3r3)  ]1//a+u>  R  -  d  t/  R  }. 
_R„  3V.nm  nm  1  '  J  nmnm'nmJ 


(5-2.7) 


n  m 


Thus  we  have 


2(-i)ntl’1  KiJl  rf(^)-i 

ni  A  9  ^tJ  ' 


(5-2.8) 


where 


f(uj)  =  [  R  oj3  -  2  d  t id  +  (t3  -  c1_3r3)]1//’*  +  uu  R  -  d  t/  R  . 

nm  nm  1  n  rr>  n  m  '  m  m 


n  m  n  m  '  n  m 

(5-2.9) 


Since  \  =  Oatu  =  u1,  X  =  0  at  oj  =  ujb  »  and  Ia  can  be  written  as 


T  2  (-1)  n+m  K,  n  .  fA,  «•  -  ndm  1 

11  -  "i  A  109  L - fW - J‘ 


(5-2.10) 


Eliminating  the  "i",  the  final  expression  for  the  real  contribution 
co  I3  is 


n+m 


2  ( -1)  “T‘"  f  n  ,  .  ,  fA,  “»  -  ndm  T/nRrr.1  l  ...... 

I.  -  In.  A  log  [  jt— ; - „  J  .  (5-2.11) 


n  m 


The  expression  for  RP  is  then  the  sum  of  Ij  ,  equation  (5-2.4), 
and  Ia ,  equation  (5-2.11),  as  follows: 


2 ( -1) n+m 
TT 

r 

Re  - 

[uai  (Kn 

yi 

2(-l)n+m 

Im 

(k,  n 

loc  r nRm 

TT  R 

lKl 

log  - 

n  m 


n  m  'n  iuTl 

T 3TJ - J/  f 


(5-2.12) 


where  u  =  x  +  iy. 
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We  have  avoided  the  singularity  of  nP^  by  setting  the  numerator 
zero  when  the  denominator  is  zero  to  obtain  a  finite  limit. 

Numerical  calculations  have  shown  that  Ij'  is  a  reasonably  well 
behaved  function  so  that  lx  can  be  evaluated  using  standard 
quadrature  methods.  Some  difficulty  is  encountered  in  the  immediate 
vicinity  of  the  singularity  at  t  =  ci  -1  R^,  but  i/  is  smooth 
enough  to  obtain  sufficiently  accurate  results  even  in  this  range. 
The  pressure  can  now  be  evaluated  numerically  with  little 
difficulty  using  complex  arithmetic. 
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APPENDIX  A 

MATHEMATICAL  DETAILS  IN  THE  INTEGRATION  FOR 

THE  STEP  WAVE  RESPONSE  P 

n  m 

Many  mathematical  details  in  the  integration  of  the  step  wave 
pressure  response  ^p  have  been  disregarded  in  the  interest  of 
arriving  at  a  solution  of  the  problem  with  as  little  complication 
as  possible.  But  many  of  the  mathematical  details  are  essential 
foe  a  good  understanding  of  the  problem  and  are  thus  the  subject 
of  this  appendix. 

A-l  CONCEPTS  AND  THEOREMS  FROM  COMPLEX  ANALYSIS 

Rosenbaum  has  used  the  "method  of  residues"  to  evaluate  the 
complex  integral  for  nPm#  This  procedure  utilizes  two  important 
theorems  from  complex  analysis:  (a)  Cauchy's  theorem  and  (b)  the 
residue  theorem.  A  few  preliminary  concepts  in  the  theory  of 
complex  variables  will  be  discussed  first. 

The  set  of  complex  numbers  is  composed  of  numbers  z  ■  x  +  iy 
where  x  and  y  are  real  numbers  and  ia  *  -1.  Complex  numbers  are 
an  extension  of  the  real  numbers  and  can  be  represented  in  a 
plane  with  x  and  y  coordinates  as  shown  in  Figure  A-l.l  below. 


FIG.  A-l.l  THE  VARIABLE  z  IN  THE  COMPLEX  PLANE 
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The  complex  number  z  can  also  be  represented  in  polar  coordinates 
(r,  9)  where  r  =  '  ul  =  (x2  +  is  the  absolute  value,  or 

modjl  is,  of  u  and  P  =  arctan  (y/x)  is  the  argument  of  u.  The  complex 
number  z  in  polar  coordinates  is  then 

z  =  r  exp(iQ)  =  r(cos  0  +  i  sin  0).  (A-l.l) 

An  important  class  of  functions  of  a  complex  variable  is  the 
"analytic  function".  A  function  f(z)  is  termed  "analytic"  in  the 
interior  of  a  closed  curve  if  that  function  is  defined  and  has  a 
derivative  at  each  point  inside  the  curve.  The  function  must  also 
be  single-valued  in  the  area  of  definition,  that  is,  for  each  z  the 
function  f(z)  mus»  ;  have  one  and  only  one  value. 

Those  points  ac  which  the  function  is  not  analytic  are  termed 
"singularities".  If  the  singularity  can  be  enclosed  by  a  circle, 
however  small,  which  does  not  include  other  singularities,  it  is 
an  "isolated"  singularity.  Singularities  such  as  those  on  the 
"branch  line"  of  a  complex  square  root  are  not  isolated. 

We  are  now  ready  to  state  the  two  theorems  used  in  the  "method 
of  residues".  Cauchy's  theorem  states  that  if  a  function  f(z)  is 
analytic  inside  and  on  a  closed  integration  path  C,  then  we  have 

f(z)  dz  =  0.  (A-1.2) 

J  C 

The  residue  theorem  states  that  the  integral 

f  ( z)  dz 
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around  the  closed  path  C,  which  encloses  only  isolated  singularities, 
is  equal  to  2ni  times  the  sum  of  the  "residues"  of  the  enclosed 
singularities.  The  counterclockwise  direction  along  the  path  C 
is  taken  as  positive.  If  the  function  f(z)  has  N  singularities 
enclosed  by  C,  then  we  obtain 

N 

f(z)  dz  =  2ni  Rj  ,  (A-1.3) 

C  j  =  l 

where  Rj  is  the  "residue"  of  f(z)  at  the  jth  singularity. 

The  "residue"  of  a  function  f(z)  at  an  isolated  singularity 
z0  can  be  defined  as  the  coefficient  aj  in  the  power  series 
expansion 

00  00 

f(z)  =  '(’z~^mz0Yn  +  150  +  Z  bn^Z  “  zo)n  •  (A-1.4) 

n=l  n=l 

This  expansion  is  called  a  "Laurent's  series  of  f(z)  about  z0". 

In  applying  the  above  theorems  for  the  evaluation  of  P  ,  we 

n  m 

must  examine  the  functions  in  the  integrand  P*  and  the  location 

n  m 

and  nature  of  their  singularities.  The  integrand  nP^  contains  square 

roots  which  are  multiple-valued  functions  in  the  complex  plane.  We 

must  then  define  single-valued  functions  in  the  region  of  integration. 

To  illustrate  the  procedure  for  defining  a  single-valued  function 

l/a 

from  a  multiple-valued  one,  we  will  study  the  square  root  z  as 
a  simple  example  of  a  bi-valued  function. 

A-2  STUDY  OF  THE  SQUARE  ROOT  w  =  zl//a 

l/2  • 

The  complex  square  root  w  =  z  '  is  defined  in  polar  coordinates  as 

zl/a  =  rX//aexp(j)  =  rl//a[cos  ^  +  i  sin  j  ]  . 

A-3 
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I  f  0  in  equation  (A-2.1)  is  increased  by  2tt,  z1//a  assumes  the 

negative  of  its  original  value.  To  obtain  all  possible  values  of  z1//a, 

the  argument  0  must  be  increased  by  4n,  that  is,  the  vector  shown  in 

Figure  A-2.1  must  make  two  full  circuits  in  the  z-plane.  Hence,  the 
l  /a 

square  root  z  takes  on  two  values  at  every  point  in  the  z-plane. 


The  concept  of  the  Riemann  surface  allows  for  a  representation  of 

x  /a 

such  multiple-valued  functions.  For  z  the  Riemann  surface 
consists  of  two  sheets,  R*  and  Ra ,  one  lying  above  the  other  and 
each  containing  a  single  value  of  z1^3  as  shown  in  Figure  A-2.2. 


A-4 


NOLTR  69-44 


When  the  point  2  makes  the  first  full  circuit,  say  from  -n  to  n,  the 
argument  0  =  0/2  has  changed  from  -tt/2  to  ‘HT/2  and  the  values  of  zl^a 
are  obtained  from  values  of  2  in  the  first  or  top  sheet  Rj  .  For  the 
second  circuit  from  0  ■=  tt/2  to  3n,  the  values  of  the  square  root  z1^'4 
are  obtained  from  values  of  2  in  the  second  or  lower  Riemann  sheet  R*  . 

The  transition  from  the  top  sheet  to  the  lower  sheet  takes  place 
at  the  "branch  cut"  or  "branch  line"  where  the  two  sheets  cross 
each  other  so  that  R*  merges  into  Ra  and  Rg  merges  into  Rj  ,  as 
illustrated  in  Figure  A-2.3  below, 

BRANCH  CUT 


FIG.  A-2.3  CROSSING  OF  THE  RIEMANN  SHEETS  AT  THE  BRANCH  LINE 
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i/a 

At  the  end  of  the  second  circuit  the  value  of  z  moves  from 
the  lower  sheet  R3  back  into  the  upper  sheet  Ri  .  By  this  way  an 
analytic,  i.e.  a  single-valued,  function  of  z1//a  is  defined  in  each 
sheet  except  for  places  where  the  sheets  cross  each  other  namely  on 
the  branch  cut. 

For  the  purpose  of  integrations,  the  path  of  integration  and 
the  branch  line  can  be  chosen  in  such  a  way  that  a  crossing  of  the 
branch  line  can  be  avoided.  The  function  then  remains  analytic.  The 
negative  real  axis  is  usually  taken  as  the  branch  cut,  even  though 
any  ray  originating  at  z  =  0  may  be  used.  The  choice  of  the  ray 
depends  on  the  particular  problem  under  consideration.  However,  the 
end  point  of  the  branch  cut  (z  =  0  in  our  case)  is  of  special  signif¬ 
icance  and  cannot  be  chosen  arbitrarily.  The  crossing  of  the  two 
Riemann  sheets  at  the  branch  cut  ends  at  the  poi  »t  z  =  0.  This  point 
is  called  the  "branch  point"  and  it  is  common  to  both  sheets.  That 

means .  that  the  '-oots  of  z  must  have  the  same  value.  In  our  case 
i  /a 

z  ■  0  is  a  common  root  and  constitutes  a  branch  point.  The 
other  branch  point  is  z 1/'a  -  •  and  it  is  seen  that  the  branch  line 
connects  these  two  branch  points. 

A- 3  BRANCH  CUTS  FOR  THE  MULTI PLE-VALUED  FUNCTIONS  IN  THE  INTEGRAND 

OF  P 
n  m 

Before  the  integral  for  the  step  wave  pressure  response 

can  be  evaluated  using  Cauchy's  theorem,  branch  cuts  must  be  made 

for  the  multiple-valued  functions  in  the  integrand.  Brarch  cuts 

must  be  made  for  a ^  =  (ua  +  c^-3)1^8  for  i  =  1,  2,  3,  4,  and 

V  *  {uara  +  [t  -  d^tcx  “a  +  ua)l/a]3}l/a.  The  function  ora  is  contained 

in  the  liquid  bottom  reflection  coefficient  K,  and  the  functions  (v9 

and  ar4  are  contained  in  the  rigid  bottom  K. 
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Other  singularities  of  the  integrand  which  must  be  considered 
are  the  zeroes  of  the  denominator  of  K  and  the  singularity  at 
infinity.  The  derivative  of  the  integrand  nP^  at  infinity  is  not 
defined.  The  nature  of  this  singularity  at  infinity  will  be  further 
examined  when  we  study  the  residues  of  the  integrand.  The  singularities 
arising  from  the  zeroes  of  the  denominator  of  K  are  treated  below. 

THE  SQUARE  ROOTS 

First  we  examine  the  square  roots  cr^.  Since  the  square  roots 
for  i  =  1,  2,  3,  and  4  are  of  the  same  form,  we  will  examine  only 
Qfi  =  (ua  +  ci~a)l/'a.  Factoring  orx  into  (u  +  icx  -1 ) 1  (u  -  ici”1)1^** 
it  is  clear  that  arx  is  the  product  of  two  functions  which  have  branch 
points  at  u  =-icx  -l  and  u  =  ici  -1  and  at  infinity.  In  the  following 
discussion  we  show  that  when  these  two  factors  are  combined  the 
branch  cut  can  be  made  from  -icx  _1  to  icx  -1  along  the  imaginary 
axis.  For  this  purpose  we  set  u  +  icx  ~l  •  rx  exp(iex  )  and 
u  -  icx -1  =  ra  exp(i0a)  as  shown  in  Figure  A-3.1.  Then  orx  becomes 

c*i  =  [rx  exp(i0j  )  ]1^a[r8  exp(i6a  )  ]1//a 
exp[j  (9X  +  6a)  ] 

[cos(6l  ■**)  +  i  sin  6a)1.  (A-3.1) 
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FIG.  A-3.1  THE  SQUARE  ROOTa]  =  (u2  +  C]  "2)  1/2 


For  the  square  root  [rx  exp(iQx  )  "I1  make  the  branch  cut  along 
the  ray  0X  =  -tt/2  from  icx -1  to  infinity  along  the  negative  imaginary 
axis,  and  for  [ra  exp(i03 )  ] 1  make  a  branch  cut  along  the  ray  83=-tt/2 
from  -icj -1  to  infinity  along  the  negative  imaginary  axis.  Hence  we 
define  the  principal  branch  (first  or  top  sheet  of  the  Riemann  plane) 
of  orj  as  follows: 

Qfj  is  defined  by  equation  (A-j.l)  and  is  restricted  to  the  domain 
-  |  s  0j  <  ,  “J*®8  <-T''r*  >0*  ra  >  0  ,  rj  +ra  >  2cx  -1  . 

That  this  branch  defines  an  analytic  function  except  along  the  cut 
from  -ici -1  to  icx  _1  will  now  be  shown. 

The  square  root  is  analytic  along  the  positive  imaginary  axis 
from  icx -1  to  infinity  since  no  branch  line  is  crossed.  Along  the 
negative  imaginary  axis  from  -icx  _1  to  infinity,  ax  is  analytic  since 
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both  [ri  exp(i0j )  ]1//a  and  [r8  exp(i08 )  ]x ^a  change  sign  from  VIj  to 

-  Vri  and  from  -Vfa  to  -  VFS  as  the  point  crosses  the  axis.  The  two 
sign  changes  cancel  leaving  crj  continuous.  Along  the  interval 

-  ici  _1  to  ici -1 ,  [rj  exp(i0j )  ]1//s  jumps  from  Vr*  to  -  Vr^  but 

[r8  exp(i0a ) l1/*  is  continuous  as  the  point  z  crosses  the  interval. 
Thus  aj  is  discontinuous  as  z  crosses  this  interval  and  not  analytic 
on  it.  The  square  root  orj  is  then  analytic  everywhere 
except  along  the  line  from  -icj  -1  to  ici  which  represents  the  branch 
cut . 

In  a  similar  manner  analytic  branches  of  a.,  i  =  2,  3,  4  can 
be  defined  with  branch  cuts  from  -ic^-1  to  ic^-1  along  the  imaginary 
axis.  Then  using  equation  (A-3.1)  and  its  cr^  counterparts  one  can 
show  that  on  the  top  sheets  of  the  Riewann  planes  for  on  the  square 
roots  have  (1)  positive  real  parts  in  the  right  half  of  the  u-plane, 
(2)  negative  real  parts  in  the  left  half  plane,  (3)  positive  imaginary 
parts  in  the  upper  half  plane  and,  (4)  negative  imaginary  parts  in 
the  lower  half  plane. 

A-4  BRANCH  POINTS  OF  FUNCTION  y  =  {uara  +  [t  -  dm<*l  ]a}1//a 

The  branch  points  of  y  are  the  roots  of  ya  =  0.  The  term  y*  can 
be  factored  into 

va  =  (t  -  a,  d  -  iur)  (t  -  eri  d  +  iur).  (A-4.1) 

m  m 

The  solutions  are  then 

u  =  R  "a[±iTr  ±  d  (t8  -  c1'aR  "8)l/a].  (A-4. 2) 

m  m  *  m 
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With  the  substitutions  cos  X  =  d  /R  ,  sin  X  -  r/R  ,  and  t  ■  t/R 

m  m  m  m  '  m  m  '  m 

the  solutions  can  be  written 

u  =  ±  i  t  sin  X  ±  (t*  -  cx‘a)l/a  cos  X  .  (A-4.3) 

m  m  m  m 

Thus  y  has  four  branch  points  compared  to  only  two  for  crx  .  But 

branch  cuts  can  be  made  for  y  as  for  crx  between  branch  points  which 

are  complex  conjugates.  Two  cuts  are  then  necessary,  but  because 

of  the  expression  (t^  -  cx-8)1//a,  different  cuts  are  required 

depending  on  whether  r  <  Ci  _1  or  t  >  cx  ~l .  These  two  cases  must 

m  ra 

be  discussed  separately. 

A- 5  BRANCH  CUTS  FOR  y  WHEN  T  <  CX  "X 

1  m 

When  t  <  cx  -1  the  four  zeroes  of  v  may  be  written 
m 


ux  =  i  [-T  sin  l  +  (cx  ~a  -  t  *)1//a  cos  4  ]  (A-5.1) 

mm  m  m 

ux  *  i  [t  sin  X  -  (cx  “*  -  t  a)1//a  cos  X  ]  (A-5.2) 

mm  m  m 

uj  =  i  [  — t  sin  X  -  (cj  ”8  -  t _a)1^a  cos  X  ]  (A-5.3) 

mm  m  m 

ui  *»  i  [t  sin  X  +  (cx -a  -  t  M)1//a  Cos  X  ]  (A-5.4) 

mm  m  m 


The  branch  cuts  are  then  made  along  the  imaginary  axis  from  ux  to  ux 

and  from  u$  to  u*  .  Since  t  >  cx -1dm,  or  t  >  cx  _1  cos  X  ,  ux  and 

mm  m 

U*  are  below  tbe  real  axis,  and  ux  and  uij  are  above  the  real  axis. 

In  order  that  y  be  positive  for  positive,  real  u,  we  define  y  to  have 
a  positive  real  part  right  of  the  cut. 
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From  equations  (A-5.3)  and  (A-5.4)  we  see  that  ug  and  u* 

are  beyond  the  integration  path  from  ux  and  ux  on  the  imaginary  axis 

and  cause  no  problem  in  the  integral  for  P  when  t  <  cx  ”1 . 

n  m  m 

Before  we  can  carry  out  the  integration  in  this  time  range 
using  real  variables,  we  must  know  whether  y  is  real  or  imaginary 
between  ux  and  ux  .  Consequently,  we  must  also  know  whether 
|  uj  |  <  cx -1  or  |uj  |  >  cx -1  for  <  cx -1  since  crx  =  (u3  +  cx”s)1//a 

is  contained  in  y. 

A-6  y  IS  REAL  AND  |  ux  |  <  cx  _1  FOR  <  cx  ~1 

For  t  <  cx -1  the  absolute  value  I  ux  I  =  t  sin  l  -  (cx  ~a  -  t3)1^3 
m  m  m  m 

cos  l  ,  which  increases  as  t  increases,  reaches  a  maximum  at 
m  m 

ss  ci'1.  The  maximum  | ux  |  =  Ci-1sin  is  less  than  or  equal  to 
ci”1.  Thus  jux  |  is  always  less  than  cx -1  for  <  cx -1 . 

In  order  to  show  that  y  is  real  for  |  ux  |  <  cx  “x  ,  make  the  change 

of  variables  u  =  ix  and  uu  =  (cx  ”a  -  xa)X/^a .  The  following  expressions 
for  the  zeroes  of  y  are  then  obtained: 

uix  =  t  cos  l  +  (cx"a  -  t  a)X//*  sin  l  (A-6.1) 

m  m  m  m 

uua  =  t  cos  t  -  (ci'a  -  t  8 ) 1  ^a  sin  l  .  (A-6. 2) 

m  m  m  m 

Thus  y  is  rerl  when  (uu  -  uui)(ui  -  uig )  s  0,  or  equivalently,  when 

m  2  or  in  £  a a  .  It  can  be  shown  that  ux  and  ux  correspond  to 

iu  =  u>i  .  Since  |x|  *  |ux  |  <  cx  _1  on  the  integration  path  for 

t  <  cx  _1  ,  it  follows  that  uu  i  uh  .  Therefore,  y  is  real  in  the 
m 

time  range  <  cx  -1  . 
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A-7  BRANCH  POINTS  OF  v  WHEN  T  >  c,  “ 

T  m  1 

When  t  >  ci  -1  the  zeroes  of  va  are  as  follows: 
m  * 


ui 


ui 


ua 


-iT 

m 

sin  j t  + 
m 

(t  a 
m 

-  cr2)i/s  cos 

m 

(A-7.1) 

iT 

m 

sin  l  + 
m 

(T3 

-  crs) i/a  cos 

l  • 

m 

(A-7. 2) 

-iT 

m 

sin  l 

m 

(T  3 

m 

-  crv/a  cos 

m 

(A-7. 3) 

iT 

m 

sin  l 

m 

(T  8 

m 

-  C!-2)1^2  COS 

l  , 
m 

(A-7. 4) 

Since  sin  l  is  positive,  uj  and  ux  are  respectively  in  the 

fourth  and  first  quadrants,  and  ua  and  U3  are  respectively  in  the 

third  and  second  quadrants.  The  branch  cuts  are  then  made  from 

uj  to  uj  and  from  u3  to  u3  .  In  order  that  the  real  part  of  y  be 

positive  on  the  integration  paths  which  are  left  of  the  cut  from 

uj  to  ^  ,  we  define  y  to  have  a  positive  real  part  left  of  the  cut 

from  ui  to  ut  and  a  negative  real  part  right  of  the  cut. 

A-8  THE  BRANCH  POINTS  OF  y  ON  THE  TOP  RIEMANN  SHEET  FOR  T  >  cx  "X 

m 

To  evaluate  the  integral  for  P  for  t  >  Cj by  the  "method 

n  m  m 

of  residues",  we  must  know  which  branch  points  of  y  are  on  the  top 
sheet  of  the  Riemann  plane,  that  is,  which  branch  points  contribute 
to  the  integral. 

From  equation  (A-4.1)  we  may  obtain  the  following  solutions  for 
t  as  a  function  of  u  when  y  is  zero: 


t  =  a,  C  +  iur 
1  m 


(A-8 . 1) 


t  =  ox  d  -  iur  , 
m 


(A-8. 2) 
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where  crj  =  (ua  +  Ci"8)1^8. 

Since  the  integration  in  u  is  to  be  carried  out  on  the  Riemann 
sheet  which  corresponds  to  positive,  real  time  t,  we  will  show  that 
the  branch  points  of  y  on  the  top  sheet  are  the  limits  of  integration 
ux  and  uj  .  This  can  be  seen  by  considering  the  following  paths 
for  u  shown  in  Figure  A-8.1s 

1.  The  counterclockwise  infinite  semicircle  ABCDEA 

2.  The  counterclockwise  infinite  semicircle  FGHIJF 

The  paths  along  the  imaginary  axis  are  allowed  to  become  arbitrarily 
close  to  the  axis  and  extend  from  -  *  to  +  ®.  The  points  J  and  E 
approach,  respectively,  -  •  and  +  •  along  the  real  axis. 

As  u  traverses  the  above  paths,  t  from  equations  (A-8.1)  and 
(A-8.2)  traverses,  respectively,  the  paths  in  Figure  A-8.2  and 
Figure  A-8.3  below.  The  primed  letters  represent  the  image  in  the 
T-plane  of  the  corresponding  unprimed  letters  in  the  u-plane. 

Between  B*  and  c'  and  between  H'  and  g',  t  is  real.  On  the  other 
parts  of  the  paths  t  is  complex.  The  shaded  region  enclosed  by 
o'c'd'e'o'  in  the  T-plane  of  Figure  A-8.2  is  the  image  of  the  fourth 
quadrant  of  the  u-plane.  Similarly,  the  shaded  region  enclosed  by 
o'e'a'b'o'  in  the  T-plane  of  Figure  A-8.3  is  the  image  of  the  first 
quadrant  of  the  u-plane.  Thus  for  real  t  >  Cx-1dm  the  zeroes  of 
y  must  be  in  the  right  half  of  the  u-plane  or  along  the  imaginary 
axis.  Consequently,  the  only  branch  points  on  the  top  sheet  of 
the  Riemann  plane  when  t  >  cj  -1Rm  or  t^  >  cj  ~l  are  the  integration 
limits  uj  and  uj  .* 

*Tue  general  method  used  here  to  study  the  branch  points  of  y  is  that  used 
by  Cagniard  [18]  pp  65,66  to  study  his  corresponding  branch  points  for 
the  problem  in  which  both  media  are  solids. 
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A-9  y  IS  REAL  FOR  >  Cx  _1  WHEN  |u|  <  Cx  ~l 

In  the  integration  of  P  for  t  >  cx  ~l  bv  the  Mmethod  of  residues" 

n  m  m 

we  must  know  the  range  over  which  y  is  real  along  the  imaginary 
axis.  Since  «x  is  imaginary  for  |u|  >  cx making  y  complex,  the 
range  |u|  *  Cj"1  offers  the  only  possibility  for  real  y. 

Therefore  consider  y  along  the  imaginary  axis  with  >  c\  -1 
and  | u |  <  Cl“l.  The  change  of  variables  u  ■  ix  in  equation  (A-4.1) 
requires  that  for  y  to  be  real  we  must  have 

ya  =  (t  -  ®id  +  xr)  (t  -  «xd  -  xr)  >  0,  (A-9.1) 

m  ro 


and  hence 


T-CNd  +  xr  >  0  and  t  -  a.d  -  xr  >  0 
m  m 

or  (A-9. 2) 

T-orxd  +  xr  <  0  and  t  -  ff.d  -  xr  <  0  . 
m  m 

For  |x|  >  0  and  |x|r  >  0  the  appropriate  condition  on  t  is  then 

t  -  crxd  -  lx|r>0  ,  (A-9. 3) 

or  equivalently, 

t  >  (cx_a  -  xa)l/*dm  +  |x|r  .  (A-9. 4) 

The  maximum  value  of  crid  +  Ixlr  for  f  x I  <  c»  _1  is  Ct  "1R_  which 

m  m 
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occurs  when  |x|  *  cx  -l  r/Rm  <  cj  “l .  Hence  y  is  real  along  the 

imaginary  axis  for  t  >  cx_lRm  or  >  d -1  and  |x|  <  Cx -1  . 

This  completes  the  study  of  the  mtltiple-velued  functions  and 

their  singularities  in  P ' .  We  now  treat  the  singularities  of 

n  m 

P#  at  the  zeroes  of  the  denominator  of  the  reflection  coefficient 
n  m 

K. 

A- 10  SINGULARITIES  IN  THE  REFLECTION  COEFFICIENT  K 

The  reflection  coefficient  K  for  liquid  and  rigid  bottoms  are 
quite  different  and  are  to  be  treated  separately.  First,  the 
reflection  coefficient  for  a  liquid  bottom  equation  (1-7.8)  is 


K 


<*l  +  baa 


(A-10.1) 


The  singularity  occurs  if  +  ba*  =  0,  that  is  if  *  -bcr»  .  But 

from  the  above  discussion  we  have  seen  that  for  a  given  value  of  u 
the  real  and  imaginary  parts  of  and  cra  have  the  same  sign  on  the 
same  sheet  of  the  Riemann  plane;  hence,  cri  =  -bcra  is  impossible  if 
oil  and  ora  are  restricted  to  the  same  Riemann  sheets.  Thus  a  liquid 
bottom  K  has  no  singularities  on  the  top  sheet  in  which  the 
integration  is  performed.  We  therefore  conclude  that  the  integrand 
of  nPm  for  a  liquid  bottom  is  analytic  everywhere  except  on  the 
branch  cuts  and  at  infinity. 

The  reflection  coefficient  for  a  rigid  bottom  is 


ffi[(2ua  +  Ca"*8)*  -  4u*0f« cu  1  -  bora  c*.  ~4 
«iC(2ua  +  c4”a)a  -  4u*ora0f4]  +  bera  c*  ~4 


(A-10.2) 


The  denominator  of  K  has  zeroes  ±ik  called  Stonley  poles  which  can 
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be  found  numerically.  It  has  been  shown  by  Spencer  [14]  that  the 
Stonley  poles  are  beyond  the  branch  points  ±ici ~x ,  ±ic* -1 ,  and 
±ic* _1  on  the  imaginary  axis  of  the  u-plane  and  also  that  tile  are 
the  only  poles  on  the  top  sheet  of  the  Riemann  plane.  Further 
references  to  the  Stonley  poles  are  made  in  Section  3-7  of  tha 
text . 

We  have  now  completed  the  discussion  of  the  singularities 

in  the  integrand  of  nPm»  The  integration  can  be  performed  if  the 

paths  are  chosen  to  avoid  the  branch  cuts.  In  the  next  section  an 

appropriate  integration  path  is  discussed  which  can  be  used  to  reduce 

the  complex  integral  to  a  real  integral. 

A-ll  THE  METHOD  OF  RESIDUES  IN  EVALUATING  INTEGRALS 

The  application  of  the  "method  of  residue"  to  evaluate  an 

integral  in  the  complex  plane  can  be  best  explained  by  an  example. 

The  problem  of  evaluating  the  rigid  bottom  integral  for  nPm 

when  t  >  Ci_1nRm  is  appropriate.  The  integration  paths  used  are 

shown  in  figure  A-ll.l  below.  The  original  integration  path  is 

shown  as  the  curve  uiAui .  Denote  the  integrand  of  P  as  P'  . 

n  m  n  m 

From  equation  (1-11.16)  the  integrand  is 


p'  =  (-l)”41"  uKn 
n  m  tti  ffj  y  * 


(A-ll.l) 


Since  u,  cr*  ,  a3  ,  or4  ,  and  Kn  are  the  same  on  path  uj  Au*  as  on  ut  Buj 
and  since  y  has  opposite  signs  on  these  paths,  we  can  write 


n^m 


■  5  £  nP^  du- 


(A-ll. 2) 
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The  and  similar  symbols  placed  on  the  integral  sign  denote 
the  path  of  integration. 

The  integration  path  shown  in  Figure  A-ll.l  above  is  composed 
of  the  following  parts:  the  branch  line  paths  Z«  and  2*  ,  the  small 
circles  £*  ane  La  around  the  Stonley  poles,  the  paths  l for 
j  =  1  -  8  ,  and  the  large  circle  I*  which  encloses  the  other  p?ths. 
This  total  path  does  not  include  any  singularities  in  its  interior 
and  hence  by  Cauchy's  theorem  the  integral  of  nP^  over  this  path  is 
zero. 


In  the  limiting  process  the  paths  ii  and  la ,  and  ,  etc., 
approach  the  same  path  and  hence  the  integrals  along  these  paths 
cancel  in  pairs.  Paths  X*  ,  E*  and  then  become  closed  circles. 

In  paths  la  and  1*  the  vertical  parts  are  allowed  to  become  arbitrarily 
close  to  the  branch  cuts  but  do  not  approach  the  same  path  as  with 
the  because  the  multj£i.6-valued  functions  change  signs  discontinuously 
as  the  branch  lines  are  crossed.  The  curved  end  paths  of  Zq  and 
reduce  to  the  points  C,  D,  A,  and  B  at  the  end  of  the  vertical 
paths.  It  can  be  shown  that  the  integrals  are  convergent  on  these 
paths  and  hence  can  be  carried  between  the  branch  points. 

The  integration  over  the  above  path  then  reduces  to  the 
following: 


£  npm  du  +  |a  P'du  +  ja  P'du  +  P  'du  +  £  P  'du  =  0.  (A-11.3) 

Hence  we  can  write 

(A-11.4) 

„pn  -  I"' nPm  du  "  T  £  p'du  -  -  p'du  +  f>  P'du  +  #*  P'du  +  f*  p'du 

UJ 
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Path  Li  is  allowed  to  have  an  infinite  radius  and  is  said  to 
"enclose"  the  singularity  at  infinity.  Similarly,  paths  £4  and 
enclose  the  Stonley  poles  ±ik.  Applying  the  residue  theorem, 
we  obtain 


P 

n  m 


lRik +  R-ik> ]  -  if- 


P' 
n  m 


du. 


(A-11.5) 


where  R  ,  R._  ,  and  R  ..  are  the  residues  of  P#  at  «°,  ik,  and  -ik 
•  ik  -ik  n  m 

respectively.  The  negative  sign  immediately  before  (R^  +  R 
appears  because  the  contours  £*  and  are  taken  in  the  negative 
(clockwise)  direction. 

Using  the  method  of  residues,  P  can  then  be  evaluated  when 

the  residues  R^,  R^#  and  R  ^  and  the  branch  line  integral  on  5* 

are  known.  For  a  liquid  bottom  the  Stonley  poles  do  not  exist  so 
the  problem  reduces  to  finding  R^  and  the  integral  or,  Ea  .  The 
residues  for  both  liquid  and  rigid  bottoms  are  derived  in  the  next 
section.  The  branch  line  integral  for  a  liquid  bottom  is  discussed 
in  Chapter  II,  and  the  branch  line  integral  for  a  rigid  bottom  is 
discussed  in  Chapter  IV.  These  integrals  are  complicated  and  have 
different  expressions  depending  on  the  relative  sizes  of  the  velocities 

Cj  and  ca  in  a  liquid  bottom  and  c*  ,  Cg  ,  and  C4  in  a  rigid  bottom. 

A-12  RESIDUE  AT  INFINITY 

To  obtain  the  residue  Rw  of  P at  infinity  the  integrand  of 
equation  (1-11.14) ,  namely 


p/  _  (-l)n4in  _ uk2 _ 

n  m  ni  Oi  Cu^r*  +  (t  -  d^x)8]1/3  ' 


(A-12 . 1) 
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must  be  expanded  in  a  Laurent's  series  about  u  =  ®.  After  dividing 
numerator  and  denominator  of  equation  (A-12.1)  by  ua  and  substituting 
z  =  1/u,  the  integrand  becomes 

(A-12.2) 


z  , x  n+m„  , 

p'(z)  =  _ z  j _ 

n  m  ni  4i+z8Cl 


_ ■_  —L*±*nn _ t—v—) 

"a)l/a  •  ra +[ rz-d  ( l+zaCl  -» )  i A ] n Va J  * 
*  m  , 


Expansion  of  each  expression  in  a  Taylor's  series  about  z  =  0  yields 


n.  [K(z=0)+  £  a  zj1n 
nP*(z)  =  ('ni  -  j ^ - » 


1+  £  bH  z3][r  +  £  c.  z3] 

j=i  3  ra  j.i  3 


(A-12.3) 


Multiplying  the  two  series  in  the  denominator  and  dividing  their 
product  into  the  series  for  Kn,  we  obtain 


P'U) 
n  m 


(-l)n4roz 


TT1 


{i-  [K(z-O)  ]"  +  V  A  z3}  , 


(A-12.4) 


or 


p' 

n  m 


(-l)n+™  rfrCnX 

TTi  l  U 


(A-12.5) 


The  residue  is  then* 

R-  *  -  {(~i&nR-  CK(-)  :"}  .  (A-12.6) 

m 

*An  explanation  of  the  negative  sign  in  the  residue  Rw  is  given  by 
E.  J.  Townsend  in  reference  [19],  page  286. 
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For  a  liquid  bottom  comparison  with  equation  (1-7.8)  shows 

that 


lim  vn 

lim  f 

-  bora-]n 

ri. 

^_b]n 

= 

U-H» 

u-  U 

+  bera 

Li 

+  b-l  * 

(A-12.7) 


where  the  limit  was  taken  directly  after  dividing  numerator  and 
denominator  by  u.  Thus,  for  a  liquid  bottom 


(-D ri  -  b~| 

ttTr  Li  +  bJ 

m 


(A-12.8) 


For  a  rigid  bottom  L'Hopital's  rule  must  be  applied.  The  result 
is  as  follows: 


tK(-nn 


lim 

u-® 


K 


n 


lim 

u-* 


rdTT  (numerator  of  K)  *n 
^  (denominator  of  K) 


1  . 


Thus  for  a  rigid  bottom  we  obtain 


(A-12.9) 


R 


CD 


n+m 


TTl  R 


m 


(A-12.10) 


A-13  STONLEY  POLE  RESIDUES  FOR  n  =  1 
It  can  be  shown  that 


lim 

z-±ik 


(z  ±  ik)n  Kn 


is  finite  whereas 
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lim 

z-*±ik 


(z  ±  ik)m  Kn 


ia  infinite  for  m  <  n.  The  Stonley  poles  are  then  poles  of  order 
n. 

For  n  =  1  the  Stonley  poles  are  simple  poles  and  may  be 
evaluated  as  follows: 


Rik  -  i-Tk  <u  -  ik)nP;-  R-ik  -  “-ik  <u  +  ik> nPm  *  (*-1-3.1) 


It  will  be  seen  later  that  R^  *  R_ik  80  that  we  need  only  evaluate 


Rik- 


The  rigid  bottom  reflection  coefficient  is 


«i[(2u*  +  c*”*)*  -  4u*crtQT4]  +  bora  C4 

Let  K  ■  A/B  where  A  and  B  are  respectively  the  numerator  and 
denominator  of  K.  Then 


R  flim  _ ttT  (-l)m'f‘1  A _  Iflim  u(u  -  ik)  ] 

ik  "  \u-ik  [ua+Cl  ~*]x/a  [uar*  +  (T-d  a, )  8 11/8  Au-ik  B  J* 

m 


(A-13.2) 


(A-13.3) 


Denote  the  first  limit  of  R^  by  Lj  and  the  second  limit  by  La 
The  limit  Iq  is  straightforward  and  becomes 


m+1 


ttT  (A-13.4) 

(c.~,-kV/‘r(c,~,-2kV+4klllc-,-ka)l/,(c.~3-k,ll/l,'l-b(c.~,-k»)l/ac.'t 

(Cj  "a-k*)1/aC-k8r*+T*-2Td  (c»  "*-k*)  x/*+d  *  (cx  “*-ka)  ]1/* 

m  m 
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Using  the  Rosenbaum  substitutions  gx  =  (ka  -  Ci-3)1^8,  g3  =  (k8  - 

g4  -  (ka  -  C4  ~3 ) 1  ^a ,  t  =  t/R  ,  and  cos  1  »  d_/R_ ,  we  have 

mm  m  m  m 


4(-i)ro+1 

ni  »  iVCTmI-(ka-c*'acosaV-2iVl  COS  'm]}1/3 


(A-13.5) 


Substituting  a  =  t  8  -  (ka  -  cx  -a  cos8  i  )  and  f  =  4t  agi  8 cos8* 

m  m  m  m 

and  removing  the  complex  square  root  from  the  denominator ,  Lx 
becomes 


-Dm+1  {*£^7 —  *a)8  -  kagag*]-  |^r}{a  + 


i/a 


4,1-11 


rri  R 


m 


gi  (a8  +  f)1 


/1 


(A-13.6) 


Looking  ahead,  Ls  \{ill  be  imaginary;  hence,  only  the  imaginary  part 

of  (a  +  iv/?)1^8  gives  a  real  contribution  to  P  .  Using  the 

n  m 

definition  of  a  complex  square  root  and  using  the  half-angle  formula 
for  a  positive  sine,  we  obtain 


m+1 


{gx  [ 


-a 


m 


(A-13.7) 

gi(a8+f)x/8 


From  equation  (A-13.3)  ,  we  can  write  the  second  limit  as 


La 


lim  _ _  u(u  -  ik) 

u-*ih  arx[{2u8  +  C4""8)8  -  4u8er3cr4]  +  bor3  c* 


(A-13.8) 


Since  ik  is  a  root  of  the  numerator  and  denominator,  L*  is  of  the 
indeterminant  form  0/0.  Applying  L'Hopital's  rule,  we  obtain 


La 


lim  f  (u)  lim  f /  (u) 
u-*ik  g  fu)  “  u“*ik  9'  (u)  * 


(A-13.9) 
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T  lim 
La  =  u-’ik 


(A-13.10) 


_ _ 2u  -  ik _ 

^-[  (2u8+c4  ~a)  a-4uaors  ^  1+aj  |8u(2ua+c4 -a) -8uoraor4 -4u3  +  £*-)}+  aV~~c4"4 

After  carrying  out  the  limit  and  substituting  gx  ,  g3  ,  and  g4  ,  we 
find 

La  =  (A-13.J1) 


i*/4 


-ka)a-kag3g4  ]-gxk{4(St-i  -ka )  +2g3 g4  +ka  (2*-  +  3j-)  ]+ 


bk 


4g3  c* 


Now  since  the  factor  k  can  be  cancelled  from  numerator  and 

denominator  of  La  leaving  both  Li  and  L3  functions  only  of  ka  ,  we 

have  R  =  R..  . 

-ik  xk 

Combining  equations  (A-13.7)  and  (A-13.11)  yields 


(A-13.12) 


lik 


+  R  =  2R . . 
-xk  ik 


_ m 

TT  i 


t  nin+1  1 k  /(aa+f)l/a  -  all/a  P 
“  n  R  gx  l  aa+f  J  1  ' 


where 


(A-13.13) 


-ka)a-kag3  g4l-bg3/(4 


c*"4) 


-k8)'-kag3g4]-gik{4(£^  -  ka) +2g3g4 +ka  (3j-  + 


.  bk 

4ga  c4  * 
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Thus  we  have  obtained  the  Stonley  pole  residue.  The  expression  for 

P  can  then  be  evaluated  using  the  "method  of  residues"  .  The  branch 
i  m 

line  integrals  and  the  complete  solutions  for  n  =  1  are  obtained  in 

Chapters  II  and  IV.  Because  of  the  complexity  of  the  Stonley  pole 

residue,  evaluation  of  P  ior  a  rigid  bottom  when  n  >  1  using  the 

n  m 

"method  of  residues"  is  impractical.  The  method  described  in 
Chapter  V  which  uses  complex  computer  arithmetic  must  then  be  used. 
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APPENDIX  B 

CHARACTERISTIC  TIMES  IN  REFLECTION  PROBLEMS 
B-l  INTRODUCTION 

Many  important  features  of  the  bottom  reflection  process  are 
described  simply  and  adequately  by  plane  wave  theory.  The  characteristic 
times  of  the  reflection  problem  are  among  these.  Snay  [20]  has 
given  a  discussion  on  the  elements  of  the  reflection  process,  and 
a  few  of  his  points  are  briefly  treated  here.  The  reader  is 
referred  to  Snay's  paper  for  further  details.  For  simplicity  our 
discussion  is  limited  to  a  liquid  bottom,  but  the  results  are  also 
applicable  to  a  rigid  bottom.  The  discussion  in  Appendix  E  on  the 
reflection  coefficient  gives  further  insight  into  the  reflection 
process . 

In  the  plane  wave  theory  two  types  of  reflection  are  distinguished: 
(a)  regular  or  undercritical  reflection  and  (b)  total  or  supercritical 
reflection.  These  types  of  reflection  depend  on  the  angle  of  incidence 
9i  and  the  angle  of  refraction  or  transmission  03  .  The  angle  of 
reflection  is  the  same  as  the  incident  angle.  The  angles  are  determined 
by  rays  perpendicular  to  the  plane  wave  fronts  as  shown  in  Figure 
B-l.l  below.  The  angles  9j  and  9»  are  related  by  Snell's  law 


sin  9i  _  Ci 
sin  9*  -  c*  * 


(B-l.l) 
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FIG.  B-l.l  PLANE  WAVE  REFLECTION  ANGLES 

Assuming  that  ca  >  Ci  ,  then  the  refracted  angle  0a  will  be 
larger  than  the  incident  angle  9*  .  When  6a  =  90°  ,  the  refracted 
ray  lies  along  the  interface  and  no  transmission  takes  place  into 
the  lower  medium.  The  corresponding  incident  angle  is  called  the 
critical  angle,  6c# 


6  =  arcsin  (ci/c»).  (B-l.S 

c 

For  incident  angles  less  than  0c,  undercr itical  reflection,  an 
incident  pressure  history  P(t)  has  a  reflection  K*P(t)  and  a 
refracted  wave  WP(t)  which  have  the  aame  shape  as  the  incident  wave, 
expressions  K  and  W  are  respectively  the  time  independent  reflection 
and  refraction  coefficients.  The  incident  energy  is  then  divided 
between  the  reflected  and  transmitted  waves. 

For  incident  angles  greater  than  0c>  supercritical  reflection, 
the  refracted  wave  travels  along  the  interface  with  the  velocity 

B-2 
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of  the  bottom.  The  refracted  energy  is  then  returned  to  the  first 

medium  so  that  there  is  no  energy  transmitted  into  the  bottom  -  hence 

the  expression  "total  reflection"  .  In  this  case  the  reflected  wave 

is  not  an  image  of  the  incident  wave  P(t)  but  is  P(t  +  T  ) ,  the 

s 

reflection  coefficient  being  1.  The  time  T  is  the  phase  shift  and 

8 

is  related  to  the  angle  of  phase  shift  <t>  described  in  Appendix  E. 

The  supercritically  refracted  wave  is  called  the  "precursor" 
since  it  arrives  before  the  main  reflected  wave.  Plane  wave  theory 
predicts  an  infinitely  long  precursor,  but  ray  theory  is  useful 
in  describing  the  precursor  path  and  arrival  time  for  a  spherical 
wave  and  is  used  in  Section  B-2  below. 

The  figures  B-1.2  and  B-1.3  below  are  examples  of  the  pressure 
history  in  which  both  the  receiver  and  the  exponentially  decaying 
source  are  in  the  first  medium.  The  second  medium  or  bottom  is  a 
non-rigid  material  for  which  ca  >  Cj  .  The  curve  of  Figure  B-1.2 
represents  undercritical  reflection,  and  the  curve  of  Figure  B-1.3 
represents  supercritical  reflection. 

In  the  regular  or  undercritical  reflection  the  direct  wave 
arrives  first  followed  by  the  reflected  wave  which  is  an  image  of 
the  direct  wave.  There  is  no  precursor. 


DIRECT  REFLECTED 

ARRIVAL  ARRIVAL 


FIG.  B-l  .2  UNDERCRITICAL  REFLECTION 
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MAIN 


FIG.  B-l  .3  SUPERCRITICAL  REFLECTION 

In  the  supercritical  case  illustrated  above  the  precursor  arrives 
first  followed  by  the  direct  wave  and  then  the  main  reflection. 

In  other  cases  the  direct  wave  might  arrive  before  the  precursor, 
but  in  any  case  the  precursor  arrives  before  the  main  reflection.  The 
■hape  of  the  reflected  wave  is  different  from  that  of  the  direct  wave. 
B-2  ARRIVAL  TIMES  OF  DIRECT  WAVE  AND  SIMPLE  REFLECTION 

In  Figure  B-2.1  below  representative  paths  are  shown  for  the 
direct  wave,  the  main  or  simple  reflection,  and  the  precursor.  From 
this  diagram  the  derivation  of  the  arrival  times  of  the  direct  wave 
and  the  main  reflection  are  trivial.  They  are  as  follows:  The 
direct  wave  arrival  time  is 


(B-2.1) 


where  R  =  (z8  +  ra)l/8. 


The  arrival  time  of  the  main  reflection  is 


i  Ti  -  ci  1  i  Ri  » 


(B-2. 2) 
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where  x  Rj  =  C(2h  -  z)a  +  ra]1,/a. 


FIG .  B-2 . 1  PATHS  OF  WAVES  WITHOUT  SURFACE  REFLECTION 
In  the  spherical  wave  problem  the  arrival  time  of  the  main 
reflected  wave  (n,m)  is  similarly 


T 

n  m 


ci _1  R 
n  m 


(B-2. 3) 


where  R  =  (  da  +  ra)1//a. 
n  m  n  m 

The  peak  of  the  main  reflection  occurs  at  the  instant  the  wave 
arrives  at  time  t  =  ci -1  R  .  The  pressure  P  at  t  s  c.'1  R 
cannot  be  obtained  directly  from  the  equations  of  Chapters  II,  IV, 
and  V,  but  must  be  obtained  from  equation  (1-11.14)  by  a  limiting 

process  as  t  -  ci -1  R  .  It  has  been  shown  (for  example  see  Towns 

n  m 

[21])  that  for  an  angle  of  incidence  6  less  than  the  critical  angle 
0C  the  peak  of  nPm  (In  the  unit,  u.ed  in  thi.  paper)  i.  given  by 
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P 

n  m 


K  /  R 
p  n  m 


at  t 


n^m 


(B-2.4) 


where  is  the  plane  wave  reflection  coefficient  given  by  equations 

(E-4)  and  (E-8)  .  But  for  angles  of  incidence  above  critical,  P 

n  m 

tends  logarithmically  to  infinity  as  t  *•  cx  R  . 

n  m 

The  arrival  time  of  the  precursor  is  more  complicated  than 
those  of  the  direct  wave  and  the  main  reflection  and  is  thus  treated 
in  more  detail  below. 

b-3  arrival  time  of  the  precursor  from  ray  theory 

Using  ray  theory  the  wave  which  we  call  the  precursor  can  be 

visualized  (see  Figure  B-2.1)  as  composed  of  rays  which  have  incident 

angles  6p  between  0c  and  9.  These  waves  travel  along  the  interface 

with  the  velocity  c*  of  the  bottom  and  finally  leave  the  bottom 

at  the  angle  6  . 

P 

For  the  case  cx  >  c»  there  is  no  real  critical  angle  and  no 
precursor.  A  wave  traveling  along  the  interface  with  the  slower 
velocity  of  the  bottom  arrives  after  the  simple  reflection. 

The  path  of  the  precursor  wave  shown  in  Figure  B-2.1  corresponds 
to  the  precursor  wave  n  =  m  =  1  of  our  spherical  wave  solution.  The 
first  response  of  the  precursor  travels  along  the  path  with  incident 
angle  0  .  The  arrival  time  6  =  x  6X  of  the  precursor  wave  for  n  =  m  =  1 
is  then 


6  =  cx  "l  (xx  +  xa  )  +  c»  _1  X3  . 


(B-3.1) 
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From  Snell' 8  law  we  have  ein  9  --  cx /c8  *  C8-1/ci-1;  hence 


cos  0c  =  (cx  ~3  -  c,"a)l/a/crl. 


(B-3.2) 


Now  xx  =h  sec  9  ,  xa  =r-h  tan  0  - (h-z)  tan  0  ,  and  Xs  =  ( h - z ;  sec  0  .  Then 


6  =  ci (2h  -  z)  sec  9  +  c8-1[r  -  (2h  -  z)  tan  9  ].  (B-3.3) 

c  c 


Substitution  of  expressions  for  sec  9  and  tan  9  yields 

c  c 


*  -  ^2h-:-  siny.  *  c* -1  [r  -  ^ 


(CX  -  -  C,  ■) 


(Ci"a  -  ca-a)V= 


6  =  r  c."1  +  (2h  -  z)  (cx  a  -  c8‘s)1/-. 


-a  _  -axi/a 


(B-3.5) 


In  this  spherical  wave  notation  the  height  from  the  image  source 
to  the  receiver  is  i  dx  =  2h  -  z  so  that  x  Rx  =  (xdx*  +  ra)1//a 
and 


i  =  r  Cj  1  +  xdx  (cx  3  -  c8  a)ll^a,  (B—  3.6) 

With  sin  =  r/x  Rx  ,  cos  =  l  dx /x  Rx  ,  and  6^  =  6/x  Rx  we  can  write 

£>m  =  c* -1  sin  tm  +  (cx“a  -  c8~a)l/a  cos  (B-3.7) 

The  above  equations  give  6  for  the  pressure  response  x  Px  .  The 
precursor  arrival  time  ^6^  for  can  similarly  be  derived  using 

ray  theory  and  can  be  written 
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6  a  r 
n  m 


+  d 


n  m 


(ci‘a 


cs-a)l/a 


(B-3.8) 


The  above  derivation  of  ^6^  uses  only  geometrical  ray  considerations. 
We  will  show  below  that  the  same  equation  is  obtained  from  equation 
(1-11.14)  for  nPm*  Hence  geometrical  ray  theory  correctly  determines 
the  arrival  times  of  not  only  the  direct  wave  and  the  main  reflection, 
but  also  the  precursor. 

B-4  ARRIVAL  TIME  OF  THE  PRECURSOR  FROM  THE  EQUATION  FOR  P 

n  m 

Aside  from  ray  theory,  the  arrival  time  ^6^  is  interpreted 
as  the  time  for  which  ^P^  first  becomes  non-zero.  Since  ray  theory 
predicts  there  is  no  precursor  if  ci  >  ca ,  we  will  first  assume  that 
c»  >  ci  and  also  that  t  <  ci  -1  R  and  c8  1  <  I  u,  I  <  ca  _1 .  The 
integration  path  is  then  that  shown  in  Figure  B-4.1  below.  The 
placing  of  ui  below  the  real  axis  and  ui  above  the  real  axis  can 
be  shown  to  be  correct  by  using  the  result  that  nP^  =  0  for 
t  <  ci  1  d  .  The  choice  of  lui  I  <  Ci -1  is  not  arbitrary.  In 
Appendix  A  it  is  shown  that  when  t  <  Ci_1nRm  then  | ui  |  <  Ci  1 . 

It  must  be  remembered  that  ±ici -1  and  ±ic3-1  are  respectively 
branch  points  of  »i  =  (u8  +  ci  a)1//a  and  aa  =  (ua  +  c3  a)1//a. 

Also  ui  and  are  branch  points  of  y  =  |uara  +  [t  -  ndm(ci~a  +  ua)^]aj^. 

When  t  <  ci  _l  R  the  square  root  (t  -  cx -1  R  )1//a  is  imaginary 
n  m  n  m 

so  that  Ui  and  ui  are  also  imaginary.  Thus  the  integration  path  is 
taken  arbitrarily  close  to  the  imaginary  axis  in  the  right  half 
plane.  The  end  points  Ui  and  Ui  are  singularities  of  the  integrand, 

>  being  zero  at  these  points.  However,  the  integral  converges  and 
the  integration  can  be  carried  from  m  to  Ui  . 
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i  lm(u) 

0) 

(2) 

_  M«j) 

13) 

(4) 

FIG.  B-4.1  INTEGRATION  PATH  IN  THE  u-PLANE  FOR 

T<  cl'VmAN0  c2  >C1 

As  the  path  approaches  the  imaginary  axis,  the  end  paths  vanish 
and  make  no  contribution  to  the  integral.  The  square  roots  crj  ,  or8  , 
and  v  ate  all  real  and  positive  on  paths  (2)  and  (3)  since  |u|  <  cj  ~l 
and  |  u |  s  ca  -1 .  The  change  in  vari*.bles  u  =  ix,  for  real  x,  in 
equation  (1-11JL4)  for  nPm  yields  for  the  contribution  I  along 
paths  (2)  and  (3) 

(B-4.1) 

I  =  -(-l)n+™  ±x  J°*  Kn[-x“ra+(T-ndm«,)>]',/‘dx. 

-ca 

where  ,  ora  ,  Kn,  and  v  are  now  real  functions  of  x* .  The  integral 

is  a  real  function  of  x2  only  and  is  then  zero  between  the  limits 

-Ca -1  and  c8  -x  .  A  non-zero  P  can  then  result  only  from  the 

contributions  of  paths  (1)  and  (4)  .  When  c*  >  ca  we  have  |uj  |  <  ca”1 

and  hence  P  =0  for  all  t  <  ci  -1  R  . 

n  m  ■*  n  m 
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From  the  above  conclusions  we  can  now  determine  6  .  For 

n  m 

Ci  >  ca  the  pressure  P  can  be  non-zero  only  if  t  i  Ci'1  R  . 

*  *  r  nm  1  1  n  m 

Thus  for  this  case  6  =  cj -1  R  and  there  is  no  precursor.  For 

n  m  n  m 

ci  <  c*  the  pressure  can  be  non-zero  only  if  lui  |  >  ca  -1 .  Thus 
n&m  is  the  time  t  for  which  Jui  |  =  Ca  .  Using  equation  (1-10.2) 
which  defines  ui  and  replacing  t  by  n6^  in  |ui  |  yields 


R 

n  m 


•  r  i 

Ln 


m 


d  (ci 
n  m  * 


-2  Ra 
n  m 


-  6 
n  m 


a  )  i  /a  j  _ 


ca 


-l 


(B-4.2) 


Solving  for  6  ,  we  get 
n  m 


6 

n  m 


*  r  ca 


-l 


d  (cx 
n  m'  1 


(B-4.3) 


Thus  the  same  precursor  arrival  time  ^6m  is  obtained  from  ray  theory 

as  from  the  spherical  wave  relection  theory. 

Above  we  made  no  mention  of  the  critical  angle  of  incidence  0c  , 

hence  we  must  determine  its  meaning  in  6  . 

n  m 

The  angle  of  incidence  e  of  the  main  reflected  ray  is  given  by 

0  »  arcsin  (r/  R  )  .  When  8  =  e  we  have  sin  8  =  Ci  /c8  and  6 
'nm  c  nm 

reduces  to  6  =  Ci  -1  R  the  arrival  time  of  the  main  reflection, 

n  m  1  n  m 

Hence,  if  6  is  less  than  8c<  the  critically  refracted  wav^  arrives 
after  the  main  reflection. 

Summarizing  our  results,  the  arrival  time  ^6^  of  the  first 
arriving  wave  can  be  expressed  by  the  following  equations: 

n6m  -  r  c,-‘  +  ndro(o,  -s  -  c,-a)lA  for  0  *  6C  (B-4.4) 
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and 

_&  ■  cj  _1  R  for  6  <  0  .  (B-4.5) 

n  m  n  m  c 

The  above  equation  is  also  valid  for  a  rigid  bottom  where  there 
is  an  additional  precursor  wave,  the  shear  precursor.  This  wave  is 
similar  to  the  above  precursor  except  it  propagates  along  the 
interface  as  a  shear  wave  with  velocity  c*  .  The  critical  angle  for 
the  shear  wave  and  the  arrival  time  of  the  shear  precursor  are  obtained 
by  replacing  ca  by  c4  in  the  appropriate  expressions. 

A  typical  pressure  history  for  an  exponentially  decaying  pulse 
measured  and  produced  in  a  liquid  layer  over  a  rigid  bottom  is 
shown  below  in  Figure  B-4.2.  The  case  shown  is  for  c*  >  c*  >  c* 
and  for  supercritical  reflection  of  both  sound  and  shear  waves, 
where  c9  denotes  the  sound  velocity  in  the  rigid  bottom. 

DIRECT 
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The  main  reflection  is  negative  for  this  case.  This  inversion 
of  the  peak  results  when  the  phase  shift  0  has  values  above  180° . 

Above,  we  have  determined  the  arrival  times  of  the  direct 
wave,  precursor  waves,  and  the  simple  reflection.  These  times  are 
the  characteristic  times  of  the  reflection  problem  and  are  important 
considerations  when  a  pressure-time  history  is  constructed. 
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APPENDIX  C 

REPLACING  THE  INTEGRATION  PATH  C  IN  THE  COMPLEX 

w 

t -PLANE  BY  A  PATH  ALONG  THE  REAL  AXIS 
In  Section  1-8  after  the  change  in  integration  variable  from 
u  to  t  it  was  stated  that  the  integration  path  can  be  taken  along 
the  real  axis.  But  the  change  of  variables  (u,  t)  does  not  lead 
directly  to  this  path,  and  it  is  the  object  of  this  appendix  to 
show  that  the  value  of  the  integral  is  unchanged  under  this  change 
of  path. 

The  change  in  variables  from  u  to  r  at  constant  w  is  defined 
in  Section  1-8  by  the  equation  (1-8.1) #  namely 

t  =  ®xndm  +  iur  cos  w,  (C-l) 

where  *  (ua  +  Cn“a)X//a.  The  integral  in  equation  (1-7.5)  becomes 

00 

nPm  -  (-l)n*"Re{f  f  yw  uo,,->Kn(fi)w  exp(-.T)dr]  dw}.  (C-2) 

°  C>'‘A 

In  this  equation  r  is  a  complex  variable,  and  the  integration 

is  carried  out  on  the  path  Cw<  the  exact  path  depending  on  w.  The 

path  Cw  is  shown  in  Figure  C-l  as  the  heavy  line  which  begins  at  the 

point  t  =  ci  ~x  d  on  the  real  axis  and  extends  to  infinity  in  the 
n  ra 

first  quadrant.*  The  letters  A  and  B  denote  points  at  infinity. 

♦For  further  discussion  on  these  integration  paths  see  Cagniard  [22]. 
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We  propose  to  use  Cauchy's  theorem  to  make  the  path  change. 
Thus  we  must  know  the  nature  and  location  of  the  singularities  of 
the  integrand.  Before  working  with  the  path  C ,  we  mu»t  then 
examine  the  partial  derivative  (du/dT)^  for  singularities. 


Re(T) 


FIG.  C-l  INTEGRATION  PATH  IN  THE  T-PLANE 


From  equation  (C-l)  we  obtain  the  partial  derivative 


<f?>  -  Lir  008  w  +  unV'0,‘  ] 

w 


-1 


(C-3) 


After  eliminating  u,  we  obtain 


(— ) 


w 


t  d  +ir  cos  w  [t3-ci  ”s  (  d3+r3cos3w)  ]1//a 
n  m  n  m 

[  d3  +  r3cos3w][Ta-ci  "3 (  +  r3cosiw)]1/3 

n  in  n  m 


(C-4) 


Ths  partial  derivative  (Su/St)  is  finite  on  C  for  all  w  except 

Vr  W 

w  ■  tt/2.  When  w  »  n/2  the  partial  derivative  approaches  infinity 
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at  t  =  &  1  d  .  But  the  integrand  remains  finite  since  at  these 
n  m 

values  of  t  and  w,  u  is  zero  and  the  limit  -1  ,  ^u(~) 

T  Cx  d  l  dt' 

n  m 

is  finite. 

It  is  very  difficult  to  integrate  along  the  path  C ;  hence  the 
shift  of  the  integration  to  the  real  axis  is  necessary.  The  shift 
can  be  performed  using  Cauchy's  theorem  if  we  close  the  path  Cw 
with  the  infinite  arc  A  B  and  the  path  B  t0  D  along  the  real  axis. 

The  square  roots  in  (^u/^t)^  can  be  made  analytic  by  making  a  branch 
cut  from  -t0  to  t0  along  the  real  axis,  where  t0  =  Ci -1  (n<^  +  r*co®aw)^« 
The  integration  path  from  t0  to  D  is  then  taken  above,  but  infinitely 
near,  the  cut.  The  integrand,  then,  has  no  singularities  within 
or  on  the  closed  path  DABt0D.  If  we  denote  the  integrand  by  f ( t) , 
then  by  Cauchy's  theorem  we  obtain 


^  f  (  t  )  d  T 
DAB  t  o  D 


f(T)dT  «  0. 


(C-5) 


The  arc  A  B  at  infinity  makes  no  contribution  since  the  integrand  is 
identically  zero  there.  Then  the  integrals  reduce  to 


~X  <3 


£w  f  (  t )  d t  =  -  J  n  m  f  (  t)  d t  =  J  f  (  t )  d t  .  (C-6) 


Ci  A  d 

1  n  m 


The  integration  for  nPm  is  now  performed  along  the  real  axis. 


and  equation  (C-2)  reduces  to 


n*m 


( _i)  n+inRe 


|2 

cn 


l"/2r 


J  -1 

o  Ci 


d 

n  m 


u<>i-1Kn(|^-)  exp(-ST)  dTdwj-.  (C-7) 
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This  is  the  desired  equation  since  it  allows  us  to  identify  the 
second  integral  as  a  one  sided  Laplace  transform. 
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APPENDIX  D 

CONVOLUTION  INTEGRAL  FOR  A  STEP  WAVE  SOURCE 
The  convolution  integral  which  allows  us  to  obtain  pressure-time 
histories  from  step  wave  solutions  can  be  derived  from  the  convolution 
integral  theorem  of  Laplace  transforms.  This  theorem  states:  If 
q,  f,  and  are  the  Laplace  transforms  with  respect  to  t  of  the 
functions  q,  f,  and  g  such  that  rj  =  f  •  g,  then  we  have 

q  =  J  f(v)g(T  -  v)  dv  *  J  f(T  -  v)g(v)dv.  (D-l) 

o  o 

The  transform  P^  of  the  pressure  response  P^  to  an  arbitrary 
source  =  f(t  -  cj-1R)/R  can  be  expressed 

P^^  =  sF(s)  A^(r  ,z,s)  ,  (D-2) 

where  A^(r,z,s)  is  a  solution  of  the  transformed  wave  equation  (1-3.4). 
The  advantage  of  writing  P^  in  this  form  can  be  seen  by  employing 
the  convolution  theorem.  The  expression  sF(s)  is  the  transform  of 
the  derivative  f'(t)  if  F(0)  =  0.  Thus  the  application  of  equation 
(D-l)  yields 

P.  =  f'(t  -  v)  A.  (r ,z,v)dv.  (D-3) 

1  v  X 
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To  find  the  meaning  of  A^(r,z,x) ,  set  F(t)  =  h(t)  then  the 
derivative  is  f'(t)  =  6(t),  where  5 ( t )  is  the  Delta  Function*  defined 
6(t)  =s  0  for  t  /  0  and 

go 

J  &(T)dT  =  1. 

0 

The  pressure  response  becomes 

pi  =  |  Mt  -  v)  Ai (r  , z , v) dv  «  Ai(r,z,T).  (D-4) 

o 

The  function  A^(r,z,T)  is  thus  the  pressure  response  to  a  source 
H(t  -  cx-1R)/R.  The  value  of  the  substitution,  equation  (D-2) ,  is 
evident.  Once  solutions  A^(r,z,T)  for  a  source  =  H(t  -  cx-1R)/R 
have  been  obtained,  then  solutions  for  a  source  ld  =  F(t  -  cx”xR)/R 
can  be  obtained  by  using  equation  (D-3)  . 

As  an  example,  consider  the  exponential  pulse 

F ( t  -  cx"xR)  *  P0H( t  -  cx"xR)  exp[  - (t  -  cx”xR)/0],  (D-5) 

where  P0  is  the  peak  pressure  at  unit  distance  from  the  origin  and 
9  is  the  time  decay  constant.  Then  we  obtain 

(D-6) 

F#(t  -  v)  «- 1®-  HIt  -  v)  exp[-(t  -  v)/e]  +  6(t  -  v)  P0  exp[-(T  -  v)/6] 


♦Good  discussions  that  make  the  delta  function  plausible  are  given 
by  Mackie  [23]  pp.  13,14  and  Churchill  [24]  pp.  26,27. 
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and 


(D-7) 


Pi(*r)  =  Po{Ai(r»z,T)  -  i  J  H(t  -  v)  exp[  - (t  -  v)/©1!  Ai(r,z,v)dv}. 


The  pressure-time  history  for  an  exponential  pulse  source  may  then 
be  obtained  from  the  step  wave  solution. 

The  step  wave  solution  in  the  water  Aj  (r,z,T)  may  be  represented 
using  the  notation  introduced  in  Sections  1-4,  1-5,  and  1-6  as  the 
following  sum: 

4  CD 

A,  U. 2.-0  -  .Pi  +  oP,  +1  l  nPm.  (D-8) 

m=l  n=l 


where  0Pi  and  0Pa  are  respectively  the  surface  reflection  and  direct 
wave  for  the  step  wave  source  and  are  derived  in  Appendix  F.  At 

any  time  t  only  those  P  for  which  t  >  6  are  non-zero.  Applying 

the  convolution  integral,  equation  (D-3) ,  to  equation  (D-8) ,  we 
obtain  the  total  pressure  Pj  (t)  in  the  water  for  a  source 
Pd  =  F (  t  -  cx  -1R)/R: 

4  oo 

pi  (T)  =  PS(T)  +  pd(T)  +  ^  I{  J  F  v*  nPm*V*dv}'  *D"9* 

m=l  n=l  0 

where  Pg  and  P^  are  respectively  the  surface  reflection  and  the  direct 

wave.  Thus  the  solution  in  the  water  for  an  arbitrary  source  P^  may 

be  obtained  from  the  sum  of  step  wave  solutions  P  . 

n  m 
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APPENDIX  E 

THE  PLANE  WAVE  REFLECTION  COEFFICIENT 
IN  THE  SPHERICAL  WAVE  SOLUTION 

The  reflection  coefficient  K  obtained  in  the  spherical  wave 

solution  can  be  identified  as  the  plane  wave  reflection  coefficient 

Kp  for  some  values  of  the  integration  variable  u.  The  rigid  bottom 

reflection  coefficient  K,  equation  (3-6.2) ,  is 


K  =  tti[(2u”  +  c*la)a  -  4u^orac»4]  -  bOf3  c« 
crx[(2u3  +  C4  a)a  -  4uaa3Qf4]  +  bo3  C4 


(E-l) 


For  a  liquid  bottom  c4  =  0,  and  ca  and  o/a  replace  c3  and  cr3  .  The 
above  equation  then  reduces  to  the  liquid  bottom  reflection  coefficient 


K 


_  on  -  btta 
«i  +  bar* 


The  variable  u  can  be  expressed 


(E-2) 


u  *  icj  -1  sin  0 , 


( E— 3) 


where  0  can  be  interpreted  as  the  angle  of  incidence  of  a  plane 
wave  for  imaginary  u  for  which  |u|  f  cx -1  .  Substitution  of  this 
expression  for  u  in  K  yields  the  plane  wave  reflection  coefficient 
Kp  which  can  be  put  in  the  following  form  given  by  Snay  [20]: 


*1  ±  Kj  -  Kj 
Kx  +  Ka  +  Ks  ' 


(E-4) 
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where 


Ki  =  cos  0  [1  -  2(c*/ci)a  sin8©]8  (E-5) 


Ka  -  4(c4/ci  ) 3  ain80  coa  0  Vl  -  (C4/C1 )  asin80  (E-6) 

*J  1  -  (c*/cj )  aain8ft 


K«  «  (PiCj/paC*)  Vl  -  (ca/ci )  aaina0  .  (E-7) 


For  a  liquid  bottom  we  have  C4  *  0  and  c*  ■  Ca  .  The  above  expresaion 

of  K  reducea  to 
P 


K  .  CO"  9  -  K* 

*p  coa  0  +  Ka  * 


(E-8) 


The  plane  wave  reflection  coefficient  Kp  given  by  equationa 

(E-4)  and  (E-8)  ia  real  with  abaolute  value  leaa  than  one  if 

0  <  9  =  arc8in  (cx/c*).  For  aupercritical  0  (0  a  0  )  the 

cr  c* 

reflection  coefficient  ia  complex  aince  the  aquare  root 

Jl  -  (ca/cj )  8ain80  i8  imaginary.  The  coefficient  Kp  can  then  be 

repreaented  in  the  complex  plane  aa  ahown  in  Figure  E-l  below. 

The  angle  0  is  called  the  "phaae  ahift"  and  ia  defined  by 

Im(K  ) 

tan  0  ■  -  —  ■  P  ■  .  (E-9) 

Re(Kp) 

The  negative  aign  in  the  phaae  ahift  ia  choaen,  when  the  aquare  roota 
are  taken  to  be  positive ,  so  that  the  plane  wave  solution  will  remain 
finite  aa  the  depth  z  below  the  source  approaches  infinity. 


NOLTR  69-44 


IMAGINARY  AXIS 


FIG.  E-l  SUPERCRITICAL  PLANE  WAVE  REFLECTION  COEFFICIENT 

For  a  liquid  bottom  when  0  *  e>cr  the  absolute  value  or  modulus 

of  K  is  one,  and  the  phase  shift  increases  from  0°  at  0  *  0 
p  cr 

to  180°  at  9  *  90°.  However,  for  a  rigid  bottom  a  second  critical 

angle,  6  =  arcsin  (ci/c*)#  which  is  greater  than  9  affects  the 

cr  s  cr 

phase  shift  and  modulus  of  K  .  For  9  =  0  _  the  coefficient  K 

P  cr  p 

has  a  modulus  of  one  and  a  zero  phase  shift.  In  the  range  9  <  0  <  0  , 

r  *  cr  crs 

|Kp|  <  1  and  the  phase  shift  takes  on  values  0  *  0*180°.  When 

9  st  e the  modulus  Ik  I  is  again  one  and  the  phase  shift  0  increases 

crs  p1  J 

to  a  maximum  value  180°  <  0  <  360°.  Depending  on  the  velocities 

max  J 

of  sound  and  shear  propagation  and  densities  of  the  two  media,  the 
phase  shift  0  can  increase  to  values  well  over  180°  in  contrast 
to  the  liquid  bottom  which  has  a  phase  shift  of  at  most  180°. 

A  simple  interpretation  of  u  in  terms  of  9  can  be  made  in  the 
precursor  integral  for  nPm  along  the  integration  path  shown  in 
Figure  2-1.1.  From  equations  (1-11.16)  and  (2-1.5)  the  change  of 
variables  u  =  ik  can  be  seen  to  yield  the  following  form  of  the 
precursor  integral: 


E-3 
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P 

n  m 


-(-1) 


n+m  2 


2_  MUl  '  x_ 

rri  J  -i  otx 
Ca 


K 


n 


dx, 


(E-10) 


where  orj  ,  y,  and  Kn  are  now  functions  of  x  and  (uj  |  =  nR~a  C^r  - 

d_  (cj  “s  R*  -  The  change  of  variables  x  =  cx  _1  sin  6 

n  in  n  m  J 

yields 


P«i 

P  =  f(S)  dl  , 
n  m  .1  « 


(E-ll) 


cr 


where  0X  =  arcsin  (ci  |ux  |).  The  angle  9X  increases  from  9cr  when 

t  =  6  ,  the  arrival  time  of  the  precursor,  to  the  "incident  angle" 
n  m 

of  the  spherical  wave,  X  =  arcsin(r/  R  ),  when  t  =  ci -1  R  . 

m  n  m  n  m 

The  spherical  wave  precursor  integral  can  then  be  interpreted 

as  follows:  The  precursor  wave  is  a  superposition  of  rays  with  angles 

of  incidence  between  9  and  0X  .  The  angle  0*  is  the  incident  angle 

for  the  slowest  ray  path  along  which  a  wave  can  arrive  at  the  receiver 

at  time  t.  As  t  approaches  ci_1nRm'  the  band  of  rays  in  the 

precursor  integral  takes  on  all  values  of  9  between  6cr  and 

The  main  reflection  (t  >  cx -1  R  )  does  not  lend  itself  to  the 

n  m 

simple  interpretation  used  above  since  the  relation  between  u  and 
9  requires  complex  values  of  0. 

In  summary,  the  reflection  coefficients  K(u)  and  K  (0)  are 

Jr 

equivalent,  and  the  spherical  wave  solution  can  be  visualized  as  a 

superposition  of  plane  waves  with  angles  of  incidence  e  which  are 

real  for  t  <  cx  1  R  and  complex  for  t  >  cx  -1  R  . 

n  m  r  1  n  m 
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APPENDIX  F 

DIRECT  WAVE  AND  SURFACE  REFLECTION 
FROM  STEP  WAVE  SOLUTIONS 


The  surface  reflection  and  the  direct  wave  for  a  step  wave  source 
are  respectively  represented  in  the  general  solution  nP^  by  the 
terms  0Pi  and  0P8  where  0di  and  od8  are  given  by  equations  (1-6.5) 
and  (1-6.6) .  Since  the  exponent  n  on  the  reflection  coefficient  K 
is  zero,  the  equation  (1-11.14)  for  0Pm  reduces  to 


”Pm-° 


T  <  6 

o  m 


Pm  "  r  u“‘"1  [uV  +  (T  -  du  T  >  n6m  • 

Ul 


(P-1) 


The  waves  0Pm  not  impact  on  the  bottom,  and  the  square  roots 

cra  or  a9  and  or4  do  not  appear  in  the  integrand.  The  branch  line 

integrals  between  -icj  _1  and  icj  -1  or  between  u*  and  uj  along  the 

imaginary  axis  are  zero.  Consequently,  0Pm  is  zero  for  t  s  cx~l0Rm. 

The  solution  0P  for  t  >  c*-10R  is  easily  obtained  using  the 

m  m 

"method  of  residues"  explained  in  Appendix  A.  Equation  (A-11.5) 

expresses  fiPm  as  the  sum  of  a  branch  line  integral,  and  contributions 

from  the  Stonley  poles  and  the  pole  at  infinity.  In  the  present 

case,  0P  ,  the  reflection  coefficient  is  K°  =  1.  Since  the  square 
m 

roots  ,  orj  ,  and  cr4  are  not  contained  in  K#  ,  the  branch  line 
integral  is  zero  and  there  are  no  Stonley  poles.  Hence,  the  only 
contribution  to  0Pm  comes  from  the  singularity  at  infinity,  and  we 
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can  write 


T 


£ 


o  P  =  -  n  i  R  T  >  Cl-1  0R 

m  ®  1  °  m 


(F-2) 


Since  K°  =  1  in  the  integrand  of  the  residue  at  infinity  R^ 

given  by  equation  (A-12.8)  or  by  equation  (A-12.10)  reduces  to 


R- 


-(-Dm 
ni  °r„  ■ 


(F-3) 


The  solution  for  0p  is  then 

m 


T  <  Cj 


m 


m 


T  >  Ci 


-I 


m 


m 


(F-4) 


The  term  0Pa  is  identified  as  the  direct  wave  and  can  be 
written 


0P3  =  PA  =  H(t  -  Cl"1  R)/R  , 


(F-5) 


where  H  is  the  Heaviside  step  function  and  R  =  0Ra  =  (za  +  . 
The  above  equation  for  PA  is  identical  to  equation  (1-4.11) .  Thus 
the  variable  t  defined  by  equation  (1-8.1)  is  indeed  the  time  t. 
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Defense  Atomic  Support  Agency 
Naval  Ordnance  Systems  Command 


This  paper  gives  a  detailed  mathematical  development  of 
a  linear,  spherical  wave  theory  for  the  reflection  of  underwater 
explosion  shock  waves  from  plane  bottoms  of  either  fluid  or  rigid 
materials.  The  Laplace  transform  method  of  L.  Cagniard  is  used  tc 
obtain  integral  solutions  for  the  pressure  which  can  be  evaluated 
numerically. 

The  paper  begins  with  linear  equations  of  motion  and  proceed) 
in  steps  through  the  derivation  of  the  wave  equations  and  finally 
to  solutions  of  the  wave  equations.  Two  methods  of  integrating 
the  integral  solution  are  discussed.  First,  the  real  part  of  the 
integral  is  separated  from  the  imaginary  part  to  allow  integratioi 
using  real  arithmetic.  Second,  a  method  of  integration  using 
complex  arithmetic  is  described.  To  aid  the  readers*  understandii  g 
of  the  problem  many  important  details  are  included  in  both 
the  text  and  the  appendices. 
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